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STUDIES  IN  PRESTRESSED 
AND  SEGMENTED  BRITTLE  STRUCTURES 

by  Ralph  L.  Barnett  and  Paul  C.  Hermann 
SUMMARY 

A  nonlinear  theory  describing  the  response  of  rectangular  pre¬ 
stressed  segmented  beams  is  verified  experimentally.  The  applica¬ 
bility  of  the  theory  is  extended  to  beam-columns  and  to  I-beams  with 
multipjle  tendons.  The  methods  of  limit  analysis  are  used  to  pre¬ 
dict  the  ultimate  load  carrying  capacity  of  prestressed  and  seg¬ 
mented  beams  and  plates.  A  statistical  hypothesis  is  developed 
and  verified  for  scaling  the  behavior  of  different  height  nonflat 
segmented  columns.  Finally,  the  feasibility  of  prestressing  cylin¬ 
drical  and  ogive  shells  is  investigated  and  evidence  is  obtained 
which  demonstrates  the  practicality  of  overwinding  as  a  prestressing 
technique . 

I.  INTRODUCTION 

To  realize  the  considerable  potential  of  ceramics  and  cermets 
in  high  performance  structures  it  is  necessary  to  circumvent  the 
problems  which  attend  brittleness  and  small  section  size.  One 
approach  to  this  problem  utilizes  the  techniques  of  prestressing 
and  segmenting,  and  indeed,  the  principal  objective  of  this  program 
is  to  study  these  techniques  for  their  possible  employment  in  aero¬ 
space  applications.  Specifically , our  goal  has  been  the  development 
of  an  analytical  capability  for  predicting  the  behavior  of  prestressed 
monolithic  and  segmented  brittle  structures  from  a  knowledge  of  the 
behavior  of  their  component  elements. 

In  the  first  phase  of  the  program,  three  fundamental  problems 
were  considered  (ref.  1)  .  The  first  of  these  dealt  with  the  devel¬ 
opment  of  transverse  tensile  stress  in  a  segmented  column  under 
axial  compressive  loading.  The  second,  involved  the  prediction  of 

the  nonlinear  response  of  a  prestressed  segmented  beam. 

- 

References  listed  at  the  end  of  text. 
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The  last  concerned  itself  with  the  benefits  which  accrue  from  pre¬ 
stressing  a  monolithic  brittle  element.  We  shall  briefly  review 
the  highlights  of  this  first  effort. 

A.  Summary  of  Previous  efforts 

1.  Transverse  cracking  phenomenon.  -  Cracking  in  a  direc¬ 
tion  transverse  to  a  uniaxial  compressive  load  was  first  recog¬ 
nized  by  F.  R.  Shanley  to  be  a  major  deterrent  to  the  application 
of  prestressing  to  segmented  members.  In  1957,  the  authors  con¬ 
ducted  a  study  of  minimum  weight  deflection  design  for  prestressed 
segmented  beams  in  which  the  roughness  of  the  segment  interfaces 
played  a  predominant  role.  Based  on  this  background,  it  was  hypoth¬ 
esized  that  the  interface  roughness  causes  transverse  cracking. 

To  support  this  view  the  following  evidence  was  established. 

(1)  The  slope  of  the  compressive  stress-strain  dia¬ 
gram  of  a  segmented  column  increases  with  in¬ 
creasing  stress.  This  is  caused  by  the  fact 
that  the  contact  area  increases  with  axial  load 
and  hence  the  stiffness  correspondingly  increases. 

(2)  Column  strength  increases  with  increasing  flatness. 

(3)  Specimens  increase  in  compressive  strength  with 
decreasing  cross  sectional  area. 

(4)  Internal  transverse  crack  lenses  can  be  observed 
in  glass  columns  (2x4xl/2  inches). 

(5)  Photoelastic  and  two-dimensional  elasticity 
results  indicate  that  an  uneven  load  distribu¬ 
tion  on  a  segment  will  cause  internal  tensile 
stresses  in  directions  parallel  to  the  interfaces. 

(6)  Triaxial  compressive  tests  indicate  a  very  sub¬ 
stantial  increase  in  axial  strength  when  a 
lateral  prestress  is  imposed. 
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2.  Load- deflection  characteristics  of  prestressed  segmented 
beams .  -  Two  quite  different  mathematical  models  were  developed 
to  describe  the  nonlinear  response  of  prestressed  segmented  beams 
with  perfectly  flat  interfaces.  The  statistical  nature  of  the 
nonflat  interface  problem  was  identified  and  its  implication  to 
both  bending  and  column  behavior  was  described.  Load-deflection 
diagrams  were  experimentally  obtained  for  segmented  glass  beams 
using  several  levels  of  prestress.  The  general  characteristics  of 
these  diagrams  are  illustrated  in  Figure  1-a  where  we  can  identify 
a  linear  and  a  nonlinear  region.  The  rough  interfaces  of  the  glass 
segments  precluded  a  deterministic  prediction  of  the  linear  portion 
of  the  curve;  however,  when  our  "perfect  interface"  models  were 
modified  to  reflect  the  proper  linear  behavior,  the  nonlinear 
region  was  predicted  with  remarkable  precision. 

3.  Strength  of  prestressed  monolithic  brittle  beams.  -  Apply¬ 
ing  Weibull's  statistical  fracture  theory,  it  was  possible  to 
theoretically  establish  for  simple  beams  a  relationship  among 
prestress  level,  reliability,  loading,  member  geometry,  and 
material  properties.  A  specific  example  was  treated  in  which  the 
prestress  results  in  a  25-fold  increase  in  ultimate  capacity  over 
a  conventional  beam  of  equal  weight  and  reliability.  The  general 
characteristics  of  the  load- deflection  diagrams  for  such  members 
are  illustrated  in  Figure  1-b  where  we  observe  the  influence  of 
both  deterministic  and  statistical  phenomena. 

B.  Summary  of  Current  Accomplishments 

The  bending  theory  previously  formulated  for  prestressed  seg¬ 
mented  beams  with  perfectly  flat  interfaces  was  verified  by  care¬ 
fully  performed  experiments  on  a  segmented  tungsten  carbide  beam. 

The  segment  interfaces  for  this  rectangular  member  were  no  more 
than  one  half  lightband  out  of  flat.  Having  placed  the  theory  on 
a  solid  foundation,  a  computer  program  was  written  to  extend  our 
analysis  capability  to  I-beams  and  box  beams  with  multiple  tendon 
arrangements.  The  applicability  of  limit  analysis  theory  for  pre¬ 
dicting  the  ultimate  load  carrying  capacity  of  prestressed  segmented 
beams  was  demonstrated  by  tests  conducted  on  a  16-foot  tubular 
alumina  beam  with  a  thin  wall  circular  cross  section. 
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Figure  1-a  Prestressed  Segmented  Bending 
Behavior 


Figure  1-b  Prestressed  Monolithic  Bending 
Behavior 
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Monolithic  Hydrostone  plaster  beams  were  used  to  accumulate  statis¬ 
tical  data  on  the  "initial  fracture"  strength  of  beam-columns. 
Weibull's  statistical  fracture  theory  provided  a  satisfactory  des¬ 
cription  of  the  measured  behavior. 

Starting  with  the  bending  theory  for  perfectly  flat  segmented 
beams,  it  was  possible  to  describe  the  behavior  of  perfect  segmented 
beam-columns  under  eccentric  axial  loads.  The  buckling  load  of  a 
perfect  segmented  column  was  shown  to  be  equal  to  the  classical 
Euler  load.  A  statistical  theory  was  proposed  for  scaling  the  com¬ 
pressive  stress-strain  diagrams  of  different  height  segmented  col¬ 
umns  with  nonflat  interfaces.  Data  obtained  for  various  size  seg¬ 
mented  glass  columns  supported  our  hypothesis  that  the  stiffness 
distribution  is  normal  and  scales  as  the  "distribution  of  the  mean". 

A  brief  investigation  into  the  behavior  of  nonprismatic  segmented 
columns  indicated  that  cracks  are  not  necessarily  arrested  at  segment 
interfaces.  Furthermore,  the  results  of  tests  on  short  plaster  back¬ 
bone  columns  suggest  that  nature  may  prefer  the  prismatic  column. 

The  theory  of  limit  analysis  was  applied  to  prestressed  seg¬ 
mented  circular  plates  and  the  resulting  predictions  agreed  closely 
with  results  obtained  from  preliminary  experiments  performed  on 
Hydrostone  plaster  disks.  Theoretically,  this  theory  provides  a 
lower  bound  to  the  strength  of  monolithic  prestressed  brittle  plates 
and  tests  conducted  on  such  elements  support  this  prediction. 

Finally,  preliminary  studies  were  conducted  with  prestressed 
cylindrical  and  ogive  shells.  The  technique  of  overwinding  was 
shown  to  be  an  effective  method  for  applying  a  prestressing  force 
over  an  extended  area.  The  first  experiments  with  segmented  plaster 
ogive  shells  seemed  to  indicate  that  the  interface  roughness  problem 
may  be  of  critical  concern. 

C.  Acknowledgments 

The  diversified  nature  of  this  program  required  the  special¬ 
ized  talents  of  many  people.  The  authors  would  like  to  acknowledge 
the  important  technical  contributions  made  by  K.  E.  Hofer  and  J.  R. 
Wingfield  in  the  area  of  experimental  mechanics,  by  J.  F.  Costello  in 
the  field  of  statistical  fracture  theory,  L.  A.  Bertram  in  the  area 
of  computer  technology,  and  by  K.  L.  Cole  in  the  area  of  filament 
winding . 
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II.  PRESTRESSED  BEAMS 

A  number  of  fundamental  investigations  involving  pre¬ 
stressed  monolithic  and  segmented  beams  are  described  in 
this  section.  Specifically,  the  theory  of  perfect  segmented 
beams  (absolutely  flat  interfaces)  presented  in  our  first 
report  is  verified  experimentally.  This  theory  is  applied  to 
prismatic  I-beams  or  box  beams  with  multiple  elastic  tendons 
and  a  computer  program  is  presented  for  establishing  the  load- 
deflection  diagrams  for  such  members.  The  possibility  of  using 
limit  analysis  methods  to  establish  the  load-carrying  capacity 
of  segmented  beams  is  briefly  exploited  with  the  aid  of  a  16- 
foot  aluminum  oxide  segmented  circular  tube.  Finally,  the  ulti¬ 
mate  capacities  of  prestressed  monolithic  hydrostone  plaster 
beams  are  measured  and  compared  to  predictions  derived  from  a 
statistical  formulation  of  the  problem  that  uses  a  beam-column 
analysis  together  with  the  Weibull  distribution  function. 

A.  Nonlinear  Bending  Theory 

1.  Theory  of  perfect  segmented  beams.  -  In  our  first  report 
two  different  mathematical  models  were  developed  to  account  for 
the  segment  separation  which  occurs  during  the  bending  of  a  seg¬ 
mented  beam  (ref.  1).  The  first  of  these,  the  incremental  model 
shown  in  Figure  2a  considers  the  beam  at  some  instant  during  the 
loading  process.  At  this  instant  the  beam  is  in  equilibrium  with 
the  applied  moment  M(x),  and  in  general,  cracks  will  have  pene¬ 
trated  into  the  beam  section  for  some  distance  along  the  segment 
interfaces.  The  relationship  between  crack  penetration  and  the 
bending  moment  at  a  station  along  the  beam  is  established  in  a 
straightforward  manner  from  moment  equilibrium.  If  an  additional 
infinitesimal  moment  6M(x)  is  added  to  this  beam,  the  resulting 
infinitesimal  response  can  be  calculated  as  the  linear  response 
of  the  uncracked  beam  section.  The  total  live  load  deflection  is 
then  found  by  summing  all  such  infinitesimal  responses  which  occur 
between  M(x)  =  0  and  M(x)  =  M(x) I  final  * 
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In  the  second  model,  the  equilibrium  model  shown  in  Figure  2b, 
the  beam  is  considered  in  its  final  loading  state.  The  portion  of 
the  beam  which  is  uncracked  is  considered  to  be  an  elastic  beam 
under  the  external  loading  M(x)  and  the  internal  loading  caused 
by  the  prestressing.  Since  the  deflection  of  an  elastic  beam  can 
be  uniquely  determined  for  every  loading,  the  deflection  of  the 
entire  beam  can  be  viewed  as  the  deflection  of  the  uncracked  portion. 

Since  the  general  demonstration  of  the  equivalence  of  these 
models  was  not  attempted  in  our  previous  study,  we  shall  deal  with 
this  problem  here. 

2.  Equivalence  of  the  two  models.  -  In  the  postulation  of  both 
the  "incremental"  and  "equilibrium"  models,  it  was  assumed  that: 

(1)  the  segment  material  is  linearly  elastic  up  to  its  ultimate 
compressive  strength,  (2)  the  interfaces  are  absolutely  flat,  (3)  the 
tendons  are  constrained  to  deflect  with  the  segments  (eliminating  any 
beam-column  action),  (4)  the  number  of  segments  is  infinite,  and 
(5)  the  resultant  prestressing  force  is  located  within  the  section 
kern  (precluding  the  existence  of  tensile  bending  stresses  and  hence, 
cracking  under  zero  external  load).  Also,  for  the  sake  of  sim¬ 
plicity,  the  equivalence  of  the  two  models  will  be  demonstrated 
for  the  case  of  rectangular  cross  section  and  zero  stiffness 
tendons . 

Both  models  are  identical  as  long  as  the  beam  is  completely 
uncracked.  In  the  cracked  region,  the  incremental  method  leads 
to  the  following  expression  for  the  deflection 


f  dAk 

A  "  Ac  +  J  “dP  dP  (  1  > 

Pc 

where  the  applied  bending  moment  distribution,  M(x),  has  been 
represented  by  Pg(x),  P  being  the  maximum  applied  bending  moment, 
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and  where  A  is  the  cracking  deflection,  P  is  the  cracking 
c  c 

moment,  and 


dAk  _ 
HP-  ' 


sa^dx  + 
s 


mg(x) 

El 


dx 


(2) 


and  where  S  and  Sc  are  respectively  uncracked  and  cracked  por¬ 
tions  of  the  beam  and  are  functions  of  P.  In  the  cracked  region, 
the  equilibrium  model  leads  to  following  expression 


A*/lrdx+/lT_dx  ( 3 ) 

S  sc  c 

where  T  is  the  resulting  moment  acting  on  a  cross  section. 

The  fact  that  these  two  models  are  indentical  will  be  demon¬ 
strated  by  showing  that 


dA  =  dAk 
clp  “dP 


(4) 


Without  loss  of  generality,  temporarily  assume  that  S(P)  is  the 
interval  0  to  a(P)  and  that  Sc(P)  is  the  interval  a(P) 
to  L. 


Then  for  Equation  (3) 


a(P) 

dA  _  _d_  f  Tm 

dP  dP  J  El 

0 


dx  4- 


d 

dP 


Tm 
,  EIC 
a(P) 


dx 


(5) 
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Thus 


,  /(P) 

dA  _  f  m_  dT  , 

Hp  J  El  Hp  dx 

0 


f° 

a(P) 


Tm 

El 


0  + 


x=0 


Tm  da 

EI  x=a(P)'^r 


m  d  /Tv, 

e  Hp  ^i  )  dx  " 

Tm 

da  . 
HP  + 

Tm 

EI 

EI 

c 

c  \ 

x=a(P) 

c 

x=L 


and  since  at  x  =  a(P) ,  I  .  =  Iy  we  have 

a(P)  L 

dA  _  f  m  dT  ,  f  m  d  /Tv, 

hp  ~  J  it  hp  dx  +  J  e  hp  dx 

0  a(P)  C 

or  more  generally 


dA  f  m  dT  ,  , 

Hp  J  ei  Hp  dx  + 
s 


/  ?  HP  (T“>  dx 

n  C 


(6  ) 


For  rectangular  cross  section  we  have 

T  bd3 
1  =  TF 


Ic  =  18  b  (I)3 


T  =  Pg(x)  -  Fe  in  S 
T  =  1  |V(e+|)  -  Pg(x) 


in  S 


(7  ) 
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Thus,  using  Equation  (7  ), 


in  S:  3p  =  g(x) 


ln  s  •  4-  (2-)  =  S(x-) 

in  bc*  3P  ''I  '  I 


Inserting  Equation  ( 8 )  into  Equation  ( 6 )  we  have  completed  the 
proof,  i.e. 


44  =  _j£  =  f  SSM 

<IP  dlF  J  El 


dx  + 


m§(x)  dx 

El 

c 


3.  General  relationships  for  multiple  tendons.  - 

a.  Arbitrary  cross  section:  The  most  general  relationships 

holding  for  any  cross  sectional  geometry  and/or  any  number'  of  elas- 

i  th 

tic  tendons,  will  now  be  described.  If  F1  is  the  force  in  the  ic 
tendon  located  e^  from  the  uncracked  neutral  axis.  Figure  3  , 
and  M  is  the  applied  moment,  then  the  resultant  moment  acting 
upon  the  cross  section  is 


where 


T  =  M  -  S  F  Q. 

i-1  x 


Qi  =  ei  +  f-  q  +  x 


and  where  n  is  the  total  number  of  tendons  - 


Define 


n  j 
F  =  Z  F 

i-1 


n  i 
2  e,- 
i-1  °  1 


Using  Equations  (11),  (12),  and  (13),  Equation  (10)  may  be  written  as 
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Figure  3  Geometry  of  Arbitrary  Cross  Section 


(14) 


n  i 

T  =  M  -  F  (e+f-  q+x)  -  2  F  (e.-e) 

i=l  1 


If  we  further  define 

n 


VL,  =  2  F1  (e.-e) 

1  i=l  1 

Then  the  expression  for  the  resultant  moment  becomes 
T  =  M-  F  (e  +  f-  q+x)-MT 


(15) 


(16) 


As  they  have  been  defined,  e  represents  the  effective  eccen¬ 
tricity  of  the  initial  total  prestress  Fq  and  M^,  represents 
the  resultant  moment  due  to  the  effective  bending  stiffness  of 


the  tendons . 


We  note  that  Mp  =  0  when  F  =  FQ  and  also  that 


Hp  =  0  for  the  case  when  all  the  tendons  have  the  same  e^  (equiva¬ 
lent  to  a  single  elastic  tendon). 

In  order  to  determine  the  crack  penetration  f  we  use  the 
fact  that  the  stress  is  equal  to  zero  at  the  top  of  the  crack. 


0  =  - 


F  4.  T  x 
A  I 


crack 

Combining  Equations  (16)  and  (17)  yields 


(17) 


M-ML, 

— ^ -  =  e  +  f  -  q  +  x  4- 


I 

Ax 


(18) 


where 

I  =  1(f),  x  =  x(f)  and  A  =  A(f). 

In  general  (except  for  the  case  of  rectangular  cross  section), 
Equation  (18)  cannot  be  solved  explicitly  for  f  .  In  the  I-beam 
computer  program  it  is  convenient  to  define  W  =  (M-Mrp)/F,  tabulate 
W  vs  f;and  use  the  table  to  accomplish  inversion  and  thus  obtain 
f  =  f(W). 

The  expressions  for  deflection  and  slope  are  simply 
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X- 

/ 


Tm 


A 

El 


TmT 


9 


9. 


I 

rl  T> 


L  J  El 
0 

L 


R  J  El 
0 


dx 


dx 


dx 


(19) 


where  the  virtual  moments  are 


rv  . 

m  =  (1  -  y^)  x  ,  for  0  <  x  <  x 


~  (L  -  x)  ,  for  x^  <  x  <  L 


(20) 


“l  =  (1  ‘  r> 


Now 


m  = 

R  L 


°o  0Lo  0Ro 


FoeL 

7eT~ 


(21) 


where  I  is  the  uncracked  moment  of  inertia, 
o 

The  equations  for  determining  the  tendon  forces  are  now,  for 
the  i^  tendon 


(Fi-F01)L 
At  Et 


Qi  <e-9o> 


x=0 


Qi  (0-eo) 


x=L 


E  J 


dx 

A 

(22) 


where  it  has  been  assumed  that  all  the  tendons  have  area  Afc  and 
modulus  Et»  There  are  n  such  relationships,  all  coupled  explicitly 
by  the  F-FQ  term  and  implicitly  by  the  other  terms,  to  be  solved 
simultaneously.  However,  due  to  the  fact  that  all  the  tendons 
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have  the  same  E^.,  these  n  equations  may  be  added  and 
solved  for  F  -  FQ  (although  transcendentally) : 


F  - 


Q(e-eo) 

„  +  Q(9-«0) 

x=0 

x=L 

rL 

L  i  n  |  dx 

A.  E.  E  J  A 
t  t  o 


(23) 


where 


(24) 


Having  used  Equation  (23)  to  determine  F  -  F  ,  the  individual 
tendon  forces  may  be  found  directly  from  Equation  (22)  . 

b.  Rectangular  cross  section:  Several  of  the  expres¬ 
sions  that  were  derived  in  the  previous  section  simplify  for 
the  special  case  of  rectangular  cross  section.  In  this  case, 
taking  dtota^  =  d  and  taking  b  to  be  the  width,  we  have  the 
following: 


A  =  b(d-f) 


T  _  b(d-f) 3 
1  IT  ■' 


(25) 


When  these  relationships  are  substituted  into  Equations  (16) 
and  (17)  we  find  that  the  uncracked  portion  of  the  beam,  S  , 
corresponds  to  the  condition  M  -  F(e+d/6)  and  the  cracked 

portion,  Sc  ,  to  the  condition  M  -  Mj,iF(e+d/6) .  Furthermore 
it  is  determined  that 


T  =  M  -  Fe  -  Mp 
A  =  A  =  bd 

O  O 

I  =  IQ  =  bd-712 


(26) 
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and  in 

Sc:  f  =  d  -  6(f) 

T  =  j  [f(e  +  j  )  +  1^,  -  m] 

A  =  6  b  (f ) 

T  3 

I  =  18  b  (f) 


(27) 


4 .  Verification  of  theory.  - 

a.  Design  of  experiments:  Of  the  various  assumptions 
entering  into  the  deflection  analysis  of  segmented  beams,  the 
most  difficult  to  realize  physically  is  that  the  segmented  inter¬ 
faces  be  perfectly  flat.  The  glass  beam  used  in  our  previous 
study  did  not  approach  this  condition  and,  consequently,  it  could 
not  be  used  to  verify  our  proposed  theory.  In  the  present  investi¬ 
gation  we  were  fortunate  to  obtain  a  set  of  100  tungsten  carbide 
gage  blocks  with  interfaces  that  were  no  more  than  one-half  light 
band  out  of  flat.  As  shown  in  Figure  4,  the  compression  load- 
deflection  diagram  for  a  18.75-inch  column  of  1-inch  x  2-inch  x 
1/4- inch  blocks  is  linear  down  to  very  low  loads  and  has  a  slope 
equal  to  that  of  a  monolithic  tungsten  carbide  bar,  that  is,  the 
modulus  of  elasticity  is  about  92  x  10^  psi. 

The  selection  of  extremely  flat  blocks  was  the  first  consid¬ 
eration  in  the  design  of  our  experimental  program.  To  minimize  the 
influence  of  the  small  range  of  nonlinearity  in  the  compression 
load- deflection  diagram  at  low  loads,  a  high  prestressing  level 
(7,000  psi)  was  chosen  for  the  tungsten  carbide  beam.  Our  previous 
work  on  glass  never  exceeded  the  1500  psi  level.  Further,  to  pre¬ 
clude  the  possibility  that  one  or  several  segments  exert  a  dispro¬ 
portionate  influence  on  the  overall  beam  behavior,  a  terminal  cou¬ 
ple-end  rotation  relationship  was  selected  for  the  response  com¬ 
parisons  . 
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Figure  4  Compression  Load-Deflection  Diagram 

for  a  18.75  inch  Segmented  Tungsten  Carbide 
Column 


17 


b.  Test  description:  The  experimental  setup  for  testing 
the  prestressed  and  segmented  WC  beam  is  illustrated  in  Figure  5. 
In  order  to  apply  terminal  couples  to  the  WC  beam,  it  was  neces¬ 
sary  to  utilize  steel  extensions  at  the  ends  of  the  beam.  Thus, 
by  employing  four  point  loading  (all  points  located  on  the  steel 
ends)  a  uniform  bending  moment  distribution  was  produced  in  the 
beam. 

The  true  (tangent)  end  rotations  were  not  measured.  Instead 
a  secant  approximation  to  the  end  rotations  was  obtained  by  meas¬ 
uring  the  relative  deflection  of  a  point  1.06  inches  from  the  as¬ 
sumed  end  of  the  beam.  In  order  to  measure  these  deflections, 
transducers  (DCDT's)  were  suspended  from  the  beam  and  located  on 
a  freely  swinging  rack,  Figure  5.  The  purpose  of  the  rack  was  to 
automatically  compensate  for  any  rigid  body  rotations  the  beam 
might  experience.  The  beam  was  loaded  such  that  it  bent  concave 
downwards  and  thus  produced  extensions  for  the  DCDT's  to  measure. 
The  experimental  setup  was  calibrated  using  a  monolithic  steel 
beam. 

Pairs  of  strain  gages  were  attached  to  each  of  the  tendons 
to  eliminate  bending  strains.  The  assembly  of  the  steel  ends  and 
the  tendons  was  calibrated  in  tension  to  verify  the  accuracy  of 
the  tendon  strain  gages.  Tendon  strain  gage  readings  were  also 
recorded  during  the  bending  tests  for  correlation  with  our  elastic 
tendon  theory. 

Two  slightly  different  beam  configurations  were  used  to  ob¬ 
tain  data.  Details  of  these  two  configurations  are  presented  in 
Figures  6  and  7.  Configuration  no.  2  is  the  better  one  from  the 
point  of  view  of  the  theory  due  to  its  simplicity.  However,  it 
has  the  significant  drawback  that  the  amount  of  preload  that  could 
be  generated  by  tightening  the  nuts  on  the  tendons  is  limited  by 
the  relatively  poor  threads  that  were  cut  on  the  tendons.  Conse¬ 
quently,  for  all  the  higher  preload  levels,  resort  had  to  be  made 
to  configuration  no.  1  which  generated  the  preload  by  extending 
the  jack  screws. 
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Figure  5  Tungsten  Carbide  Beam  Experimental  Setup 


Figure  7  Details  of  Beam  End  Configurations 


Preliminary  bending  tests  were  conducted  using  configuration 
no.  1.  Preload  levels  of  5,  8,  11  and  14  kips  were  selected  which 
were  in  the  monolithic  stiffness  range  according  to  Figure  4  for 
the  column  test.  The  tests  were  run  in  the  sequence  14,  11,  8,  5, 

5,  8,  11  and  14  without  ever  bringing  the  preload  to  zero.  The 
results  of  these  tests  are  presented  in  Figure  8. 

Two  final  terminal  couple- end  rotation  diagrams  were  deter¬ 
mined  for  the  WC  beam  where  exceptional  attention  was  devoted  to 
the  testing  details.  In  each  test,  the  preload  was  gradually  in¬ 
creased  from  zero  to  the  desired  level  and  many  load  precycles 
were  applied  to  the  beam  to  completely  stabilize  the  system.  The 
first  test  was  conducted  with  configuration  no.  1  and  a  14  kip 
preload.  The  preload  in  the  second  test  was  4  kips  which  enabled 
us  to  use  the  configuration  no.  2  and  avoid  the  complications  in¬ 
troduced  by  the  jack  screws.  The  results  of  these  tests  are  pre¬ 
sented  in  Figures  9  and  10. 

c.  Comparison  of  theory  and  experiment:  The  theory, 
with  which  the  experiments  will  be  compared,  is  in  the  form  of  a 
computer  program  which  was  developed  under  the  first  phase  of  this 
contract  (ref.  1).  In  order  to  take  into  account  the  effect  of 
multiple  elastic  tendons  and  also  to  attempt  to  compensate  for  the 
structural  complications  (especially  in  configuration  no.  1)  that 
were  introduced  between  the  end  of  the  WC  beam  and  tendons,  the 
original  computer  program  was  significantly  modified. 

Recalling  the  elastic  tendon  equation  for  the  case  of  multiple 
tendons ,  we  have 


(F1  -  Fj)  L 


Qi(9-eo) 


x=0 


Qi  <e-eo> 


F-F 


o  fk  dx 


E 


x=L 


/ 

0 


(22) 


L  Otar 

where  9  =  /  . . w  +  ■  dx 

0 


E  I 


(28) 


and  where  the  interval  (0,L)  represents  the  assumed  length  of 
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Total  Preload  In  Upper  And  Lower  Tendons  (kips) 


Figure  10  Tendon  Stress  -  Terminal  Couple  Diagram 
for  a  Prestressed  Segmented  WC  Beam 
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the  WC  beam.  In  order  to  account  for  the  steel  ends  and  (in 
configuration  no.  1)  the  jack  screws  in  Equation  (22),  the  end 
rotation  should  be  expressed  as 


L  Tim 

8  '  K1  {,  ET-  dx  +  K2  <M  -  Fe 

F  -  F 

and  the  beam  compression  term,  — g — 
fied  to 


-  V  +  K3  (  f")  (29) 

—  4^  ,  should  be  modi- 
0  A 


(F 


+ 


(30) 


where 


K, 


Assumed  Length 

1  ( 

<nr)jS 


Kc 


K, 


1  (  L  -i 

~T~  ■ ’ 


ST 


(JL1) 

''AE  ; 


JS 


Kc 


( 4-) 


ST 


(31) 


These  correction  terms  and  factors  were  derived  with  the 
following  assumptions:  (1)  M(x)  =  constant;  (2)  the  steel  ends 
have  the  same  width  and  depth  as  the  WC  beam  and  are  infinitely 
segmented  so  that  their  area  and  moment  of  inertia  depend  upon 
the  loading  in  the  same  manner  assumed  for  the  WC  beam;  and  (3) 
that  the  jack  screws  act  simply  as  two  force  members.  The  sub¬ 
scripts  ST  and  JS  refer  to  the  steel  ends  and  the  jack  screws 
respectively,  2w  is  the  vertical  distance  between  jack  screws  and 
^°JS  t*ie  extended  length  of  the  jack  screws. 
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TABLE  I 


LENGTHS  OF  JACK  SCREW  CORRESPONDING  TO  THE 
VARIOUS  PRESTRESS  LEVELS 


F0(kips) 

i°s  (inches) 

14 

0.0559 

11 

0.0439 

8 

0.0319 

5 

0.0200 

In  configuration  no.  1,  Lg,j  =  2.64  inch,  w  =  .75  inch, 

AjS  =  .0227  sq.  in.,  and  Lwc  =  20.0  inch.  In  configuration 
no.  2  LCT  =  1.44  inch  and  L  =  21.0  inch.  Also,  in  both 
configurations  Eg,p  =  Ejg  =  30  x  lO-"5  ksi,  and  the  assumed 
length  =  18.75  inch.  Insertion  of  these  values  into  Equations 
(31)  yields  Table  II. 


TABLE  II 

VARIOUS  CONSTANTS  USED  IN  THE  ROTATION  COMPUTATIONS 

FOR  THE  WC  BEAM 


Configuration  K-^ 

No. 

K2  6 
xlO-6 

K3 

xlO' 

K4 

-6  xicr6 

K5  6 

xlO  6 

14 

1 

20/18.75 

36 

44 

21 

88 

11 

1 

20/18.75 

29 

44 

16 

88 

8 

1 

20/18.75 

21 

44 

12 

88 

5 

1 

20/18.75 

13 

44 

7 

88 

4 

2 

21/18.75 

24 

0 

48 

As  mentioned 

previously,  the 

experiments 

yielded 

average 

(of  both  ends)  secant  end  rotations  corresponding  to  vertical 
deflection  measurements  taken  at  stations  1.06  inches  apart 
at  both  ends  of  the  beam.  Thus,  the  computer  program  was  ad¬ 
justed  to  compute  the  same  secant  end  rotation. 

Figure  8  illustrates  the  comparison  of  the  theory  and  the 
preliminary  experiments.  It  is  observed  that  generally,  the 
higher  the  preload,  the  better  the  agreement.  There  are  a  number 
of  reasons  why  this  is  expected;  higher  preloads  provide  a 
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tighter  system  with  greater  contact  area  between  segments, 
the  influence  of  the  low  stress  nonlinearity  is  suppressed, 
and  the  relative  errors  in  preload  determination  are  reduced. 

It  should  be  noted  that  a  rough  analysis  of  the  jack  screws 
indicated  that  they  tend  to  open  up  in  the  nonlinear  ranges 
for  the  cases  FQ  =  5  and  8  kips. 

The  comparison  of  the  theory  and  the  final  experiments 
is  illustrated  in  Figures  9  and  10.  The  agreement  is  much 
better  than  in  the  case  of  the  preliminary  experiments,  both 
in  the  linear  and  in  the  nonlinear  ranges .  This  is  due  to  the 
extra  care  used  in  performing  the  final  experiments  and  pos¬ 
sibly  to  the  fact  that  the  jack  screws  were  eliminated  in  the 
case  Fq  =  4  kips.  The  jack  screws  behaved  as  predicted  in  the 
14  kip  case. 

5.  Segmented  beams  with  nonflat  interfaces.  -  When  the 
segment  interfaces  are  not  flat,  the  contact  area  between  any 
two  segments  may  vary  from  almost  full  contact  to  almost  no 
contact.  Furthermore,  the  actual  contact  area  cannot  be  pre¬ 
dicted  for  particular  segments  since  it  varies  randomly  from 
interface  to  interface.  However,  when  all  the  interfaces  have 
been  drawn  from  the  same  population,  which  usually  happens  when 
the  same  manufacturing  technique  is  used  for  all  the  blocks,  we 
can  predict  the  behavior  of  groups  of  segments  in  a  statistical 
sense.  The  compression  test,  for  example,  furnishes  a  measure 
of  the  average  or  effective  contact  area  at  every  level  of  com¬ 
pression.  For  a  given  segmented  column  the  effective  area  is 
computed  by  multiplying  its  nominal  area  by  the  ratio  of  its 
tangent  stiffness  to  the  equivalent  monolithic  stiffness. 

The  statistical  nature  of  the  interface  contact  gives  rise 
to  a  number  of  important  implications.  First,  a  description  of 
the  compressive  stress-strain  diagram  will  be  statistical  and 
will  depend  on  both  the  column  area  and  the  length.  This  prob¬ 
lem  will  be  studied  more  thoroughly  in  Section  III.  Second, 
the  moment  of  inertia,  like  the  area,  will  be  random  with  a 
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lower  bound  equal  to  zero  and  an  upper  bound  equal  to  the  full 
moment  of  inertia.  It  follows  that  beam  behavior  will  be  statis¬ 
tical  in  nature  and  that  the  variability  experienced  will  depend 
on  the  length  and  area  of  the  beam,  the  number  of  segments,  and 
on  the  loading. 

Because  the  effective  stress-strain  curve  of  a  segmented 
column  is  area  and  length  dependent,  it  is  clear  that  an  equivalent 
nonlinear  material  cannot  be  defined  which  is  useful  for  response 
calculations.  For  example,  in  a  segmented  beam  under  sufficiently 
high  loads,  the  nominal  "uncracked"  area  diminishes  continually 
as  the  load  is  increased.  Consequently,  the  stiffness  of  the  section 
tends  to  increase  as  a  result  of  increasing  axial  stress  and  de¬ 
creasing  nominal  area;  it  tends  to  decrease  because  of  a  reduced 
moment  of  inertia.  The  theory  of  perfect  segmented  beams  accounts 
only  for  variations  in  the  moment  of  inertia.  An  exact  solution 

for  the  response  will  require,  in  addition,  that  we  describe  the 
bending  stiffness  of  any  "uncracked"  area  subjected  to  a  linear 
strain  field  that  varies  from  zero  at  one  end  to  some  arbitrary 
maximum  value  at  the  other.  It  does  not  appear  that  an  exact  bend¬ 
ing  stiffness  relationship  can  be  derived  from  column  tests.  The 
desired  relationship  can  be  obtained  from  terminal  couple-end 
rotation  tests;  however,  the  authors  doubt  whether  the  end  result 
justifies  the  effort  involved. 

An  approximate  description  of  a  load-deflection  diagram  can 
be  obtained  through  a  slight  modification  of  the  perfect  interface 
theory.  If  the  modulus  of  elasticity  of  the  segment  material  enter¬ 
ing  into  the  theory  is  replaced  by  the  tangent  modulus  of  the  beam 
associated  with  the  initial  prestress  level,  the  perfect  interface 
theory  will  predict  the  average  load  deflection  diagram.  In  our 
first  report  we  demonstrated  that  this  procedure  predicted  the  non¬ 
linear  behavior  of  segmented  beams  when  the  initial  straight  line 
portion  of  the  curve  was  matched  to  the  data.  The  necessity  of 
matching  up  the  initial  slope  rather  than  inferring  it  from  the 
tangent  modulus  taken  from  the  associated  column  tests  strikes  at 
the  core  of  the  statistical  problem.  The  tangent  modulus  prediction 
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describes  only  the  average  behavior  of  a  group  of  segmented  beams 
and  not  the  response  of  a  single  member. 

To  use  the  suggested  approximate  technique  for  predicting  beam 
response,  we  require  a  description  of  the  compressive  stress-strain 
curve  for  the  nominal  cross-sectional  area  of  the  beam.  There  is 
every  reason  to  believe  that  the  statistical  description  of  a 
single  typical  size  column  will  enable  one  to  describe  any  other 
column  of  different  area  and  length.  This  problem  is  studied  in 
Section  III. 

As  a  final  observation  concerning  nonflat  interfaces,  we  re¬ 
call  from  our  last  report  that  the  initial  portion  of  the  central 
load-central  deflection  diagrams  for  the  glass  beam  were  all 
straight.  On  the  other  hand,  the  compression  stress-strain  diagram 
for  the  glass  is  curvilinear  with  a  monotonically  increasing  slope. 
This  apparent  anomaly  can  be  explained  by  considering  an  axial 
prestressed  beam  that  is  free  from  lateral  loads.  The  initial 
bending  stiffness  of  this  member  is  proportional  to  the  tangent 
modulus  of  a  compression  curve  at  the  given  prestress  level.  Now, 
when  a  bending  moment  is  applied  to  this  member  the  compression 
fibers  will  tend  to  get  stiffer  and  the  tension  fibers  will  tend 
to  become  more  flexible.  The  two  effects  neutralize  each  other  and 
mitigate  the  influence  of  the  compression  nonlinearity. 


B.  I-Beam  or  Box  Beam  with  Multiple  Tendons 

The  analysis  of  I  (or  equivalently  box)  beams  proceeds  exactly 
as  outlined  in  the  section  on  general  relationships.  Referring  to 
Figure  11,  we  obtain  the  following  cross  sectional  properties: 

Uncracked  section  properties  (f  s  0) 
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Cracked  section  properties 
Case  I  (0  <  f  <  tf) 

A  =  b(t£-f)  +  dtw  +  bt£ 

x  =  i  [|  (t£-f)2  +  dtw(tf+|-  f)  4-  btf(|  tf  +  d  -  f)]  (33) 

t  2 

I  =  [b(tf-f)3  +  twd3  +  btf^]  +  b(t£-f)(x  -  j  ) 

2  2 
+  dtw  (x  -  ^  -  t£  +  f)  +  bt£  (j  t£  +  d  -  x  -  f) 

Case  II  (tf  £  f  <  tf  +  d) 

A  -  (tf  +  d  -  f)  tw  +  btf 

x  =  \  (tf+  d-f)2  +  btf  (|tf  +  d-f)]  (34) 

3  t 

1  =  [tw  (tf  +  d'f)  +  btf3]  +  Cw  (tf+  ”  T  -  I  +  l> 

3  -  2 

+  btf  (jtf  +  d  -  x  -  f ) 


Case  III  (t£  +  d  <  f  <  2t£  +  d) 
A  =  b (2tf  +  d  -  f) 

-  .  d  f 

x  "  tf  +  7  -  7 

K  3 

I  =  Y1  +  d  -  f) 


(35) 


Due  to  the  algebraic  complexity  of  the  above  expressions  for 
A(f ) ,  x ( f )  and  1(f),  the  equation  for  the  crack  penetration 

W=e+f-q+x+  <36) 

where 

T7  M-Mt 
"'T1 
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cannot  in  general  be  inverted  to  find  f=f(W).  In  order  to  be 
able  to  handle  this  new  complication,  a  new  computer  program 
has  been  developed.  A  listing  of  the  program  (in  Fortran  II  for 
the  IBM  7094)  plus  sample  input-output  is  presented  in  Appendix  A. 
This  program  is  currently  set  up  to  handle  I-beams  with  multiple 
elastic  tendons.  Each  of  the  cross  section  properties  A(f),  x(f) 
and  1(f)  are  programmed  as  function  subprograms  and,  consequently 
any  other  cross  section  geometry  may  be  investigated  by  merely 
changing  these  three  subprograms . 

The  inversion  of  Equation  (32)  is  accomplished  in  the 
computer  program  through  the  vehicle  of  a  table  of  W  vs  f.  This 
table  is  constructed  in  increments  according  to  input  specifi¬ 
cations.  The  inversion  is  easily  effected  by  a  function  subprogram 
which  merely  searches  the  table  using  W  and  linearly  interpolates 
between  the  bracketing  values. 

Prestressed  segmented  glass  beams  with  the  cross  sectional 

geometries  shown  in  Figure  12  and  with  E  =  10.5x10°  psi, 

6 

Et  =  30  xlO  psi,  and  L  =  Lt  =  38  inches  have  been  selected  to 
illustrate  and  compare  the  behavior  of  I-beams,  webless  I-beams, 
and  rectangular  beams  of  the  same  overall  dimensions.  Figure  13 
illustrates  the  behavior  of  these  beams  under  terminal  couples 
for  the  case  of  zero  stiffness  tendons  (Et=0) .  We  observe  that  the 
rectangular  beam  is  the  stiffest  (least  deflection)  member  in  the 
initial  load  range,  the  I-beam  in  the  middle  range,  and  finally 
the  webless  I-beam  in  the  final  range. 

It  is  interesting  to  note  that,  for  a  load  of  18  inch  kips, 
the  rectangular  beam  has  a  crack  penetration  of  2.5  inches,  the 
I-beam  has  1.7  inches,  and  the  webless  I-beam  is  still  uncracked. 
Thus,  the  web less  I-beam,  which  has  only  25  percent  of  the  area 
and  weight  of  the  rectangular  beam,  is  unquestionably  the  most 
efficient  cross  section  for  prestressed  beams  in  terms  of  stiff¬ 
ness  per  unit  weight. 
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Figure  12  Beam  and  Tendon  Geometries 
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with  Zero  Stiffness  Tendons 


We  also  observe  for  the  case  of  the  webless  I-beam  that 
although  abrupt  changes  in  character  occur  in  the  terminal 
couple  vs .  end  rotation  curve  and  the  terminal  couple  vs.  crack 
penetration  curve  as  the  crack  passes  through  the  bottom  flange, 
that  there  is  in  fact  no  discontinuity  in  slope  in  either  of 
these  curves.  One  additional  noteworthy  item  for  the  case  of 
zero  stiffness  tendons  is  the  fact  that,  as  the  cracks  penetrate 
the  cross  sections,  eventually  (f  >  +  d)  the  uncracked  sections 

become  identical  and  hence  the  behavior  of  the  three  beams  coin¬ 
cide  in  this  range. 

The  behavior  of  the  three  beams  for  the  case  of  elastic 
tendons  is  shown  in  Figure  14  and  15.  In  general,  the  effect 
of  the  tendon  elasticity  is  to  increase  the  stiffness  and  the 
slope  of  the  terminal  couple  vs .  end  rotation  curve  and  to  limit 
its  minimum  slope  as  the  load  increases  to  infinity.  There  are 
no  discontinuities  in  the  slope  of  the  curves  for  end  rotation, 
prestress  force,  prestress  moment,  or  crack  penetration  as  the 
crack  passes  through  the  bottom  flange  for  the  case  of  the  web¬ 
less  I-beam. 


As  in  the  case  of  terminal  couples  on  a  rectangular  beam 
with  a  single  elastic  tendon  (ref.  1)  it  is  expected  that  the 
cracks  will  never  pass  completely  through  the  section.  This  fact 
is  verified  by  the  crack  penetration  curves  in  Figure  15  which 
seem  to  be  approaching  asymptotic  values  other  than  ^total=  ^ ' 
Whether  all  three  curves  will  have  the  same  asymptote  depends 
upon  the  relative  stiffnesses  of  the  beams  to  the  tendons.  Thus, 
in  general,  it  cannot  be  expected  that  the  behavior  of  the  rec¬ 
tangular  beam,  I-beam,  and  webless  I-beam  will  become  coincident 
as  the  applied  load  becomes  unbounded. 
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Figure  14  Load-Deflection  Diagrams  and  Crack  Penetration  Diagrams  for  I-Beams 

T.T-I  4—  T-«  T71  1  n  n  4-  ■!  _  m  _ _ J _ _  ° 


Tendon  Moment,  Mf  (in -kips)  Total  Tendon  Force,  F  (kips) 


Terminal  Couple  (in-kips) 


Figure  15  Tendon  Force  and  Moment  vs .  Terminal  Couple 


38 


C.  Limit  Analysis 

The  possibility  of  approximating  the  load-deflection  diagram 
of  a  prestressed  segmented  beam  by  one  which  is  elastic-perfectly 
plastic  was  suggested  in  our  previous  report  (ref.  1).  Indeed, 
this  possibility  was  exploited  in  1952  by  Kooharian  in  his  study 
of  segmented  concrete  arches  (ref.  2).  In  these  arches,  the  com¬ 
pressive  forces  acting  normal  to  the  segment  interfaces  were  not 
provided  by  prestressed  tendons,  but  rather  by  the  arches'  reaction 
to  live  and  dead  loading. 

The  applicability  of  limit  analysis  for  predicting  the 
ultimate  load-carrying  capacity  of  prestressed  segmented  beams 
is  investigated  in  this  section  with  the  aid  of  the  16- foot 
segmented  aluminum  oxide  beam  shown  in  Figure  16  .  Each  segment 
in  this  member  was  four-inches  in  both  length  and  outside  diameter 
with  a  wall  thickness  of  3/16-inch.  The  beam  was  composed 
of  48  such  segments  and  the  prestressing  was  accomplished  by 
pretensioning  a  1/4-inch  steel  prestressing  tendon  that  passed 
along  the  axes  of  the  cylinders  and  was  secured  to  steel  end 
plates.  To  preclude  the  presence  of  secondary  bending  effects 
(beam- column  behavior),  the  tendon  was  constrained  to  the 
centroid  of  the  sections  by  five  closely  fitting  wooden  spacer 
inserts  on  three-foot  centers.  Two  strain  gages  on  opposite 
sides  of  the  tendon  were  used  to  monitor  the  prestress  level. 
Simple  end  supports  were  provided  by  two  saw  horses  as  shown 
in  Figure  17a  .  The  beam  was  loaded  with  dead  weights  and  the 
finest  load  increment  was  five  pounds. 

The  beam  was  tested  under  the  three  types  of  loading 
shown  in  Figure  18  .  In  all  cases,  the  loading  was  continually 

increased  until  a  0.005-inch  thick  feeler  gage  could  be  inserted 
between  the  separated  segments  to  a  depth  of  two- inches.  The 
load  associated  with  this  condition  is  recorded  in  Figure  18 
as  ^Measured*  ^he  f^rst  test  1X111  was  the  central  loading  con¬ 
figuration.  It  was  noted  at  the  conclusion  of  this  test  that 
longitudinal  cracks  had  appeared  on  the  compressive  side  of 
several  segments  that  were  close  to  the  beam's  center.  For  the 
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Figure  16  Alumina  Hollow  Circular  Prestressed  Segmented  Beam 


Figure  17  Segmented  Alumina 
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d)  Single  Quarter  Point  Loading 


Figure  18  Limit  Analysis  of  Prestressed  Segmented  Beam 
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remaining  two  tests  the  beam  was  rotated  such  that  the  final 
bearing  area  was  away  from  the  cracks.  During  the  final  loading, 
two  concentrated  loads  at  the  quarter  points,  horizontal  cracks 
developed  under  the  loads  which  resulted  in  the  catastrophic 
failure  shown  in  Figure  17b.  The  failure  modes  shown  in 
Figure  19  are  typical  of  the  primary  segment  fractures. 

The  prediction  of  the  ultimate  loading  of  the  alumina 
beam  follows  precisely  the  methods  of  limit  analysis  for  simple 
beams.  Here,  we  take  the  plastic  moment  to  be  the  prestressing 
force  times  half  the  beam  depth  as  shown  in  Figure  18a  .  The 
formulas  and  predictions  for  the  limit  loads  are  given  in 
Figure  18b,  c,  and  d  where  we  reflect  the  following  physical  data 

(1)  Prestress  level,  F  =  5000  lb 

(2)  Weight  density  of  tendon  -  0.167  lb/ft 

(3)  Span  length,  L  =  182.5  in. 

(4)  Weight  of  entire  beam,  W  =  90  lb 

(5)  Weight  of  loading  fixture  -  10  lb 

(6)  Plastic  moment,  Mp  «  (2)  (5000)  =  10,000  in. -lb 

We  observe  from  this  figure  that  the  predicted  loads  are  from 
2.92  percent  to  5.41  percent  lower  than  the  measured  loads. 

D.  Prestressed  Monolithic  Beams 

Our  previous  work  (ref.  1)  was  directed  toward  applica¬ 
tion  of  Weibull's  statistical  fracture  theory  to  monolithic 
prestressed  brittle  beams,  neglecting  the  effect  of  beam-column 
action.  That  effort  was  devoted  to  developing  the  relationships 
among  prestress  level,  load,  geometry  and  the  reliability  of  a 
structural  member.  We  then  proposed  to  demonstrate  those  results 
by  means  of  an  experimental  program  using  Hydrostone  plaster 
beams.  In  this  year's  effort,  the  additional  moment  due  to  the 
eccentricity  of  a  straight  tendon  when  the  beam  is  laterally 
loaded  has  been  considered.  This  may  well  not  be  of  major 
significance  in  a  real-life  application  where  the  prestressing 
tendon  should  be  constrained  to  deflect  with  the  beam,  but  it 
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Figure  19  Typical  Failure  Modes  of  Alumina  Beam  Segment 


is  necessary  explain  the  experimental  results  reported  in  this 
section.  The  experiments  which  were  carried  out  involved  four- 
point  bending  tests  on  the  hydrostone  beams  described  in  Appendix 
B.  The  dimensions  are  shown  in  Figure  20  ,  and  the  loading  fixture 
in  Figure  21  .  Since  all  specimens  failed  between  the  loads,  the 
experiment  can  be  considered  as  a  pure  bending  test  with  a  gage 
length  L. 

According  to  Weibull's  theory,  the  failure  probability  F 
will  be: 


F(o)  =  1  - 
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where  a  is  an  intensity  level;  o  <{  is  the  actual  stress  distribu¬ 
tion  in  the  body;  4-^,  and  £3  are  space  coordinates;  v  is  a 

unit  volume;  and  m,  a  ,  and  gq  are  statistical  distribution  para¬ 
meters.  For  convenience,  the  notation 

F  =  1  -  exp  [-B]  (38) 


is  used,  where  the  definite  integral  B  is  called  the  risk  of 
rupture. 

For  the  case  of  a  rectangular  beam  subjected  to  a  pure 
couple, the  risk  of  rupture  is 
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Figure  20  Schematic  of  Four-Point  Bending  Test 
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where  V  is  the  volume  of  a  beam  of  length  L,  width  b,  and  depth  d; 
the  maximum  fiber  stress  a  =  6M/bd2  is  taken  as  the  intensity 
level;  &  is  taken  as  2y/d;  and  y  is  the  coordinate  through  the 
beam  depth  measured  from  the  neutral  axis.  For  materials  like 
Hydrostone  that  have  low  tensile  strengths,  it  is  common  practice 
to  attempt  to  fit  their  fracture  data  with  a  two  parameter 
assumption,  i.e.,  a u  =0.  Then,  Equation  (39)  simplifies  to 
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where  the  volume  has  been  incorporated  into  the  scale  parameter 


The  cumulative  probability  of  failure  curve  for  the  unpre- 

stressed  beams  is  shown  in  Figure  22  where  the  data  have  been 

ordered  and  the  probability  of  failure  at  the  stress  associated 
th 

with  the  i  observation  is  estimated  to  be  F  =  i/N+1  where  N 
is  the  total  number  of  observations.  The  maximum  likelihood 
estimates  of  the  parameters  o'  and  m  are  found  from  the  solutions 
of  (ref.  3  ): 
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where  the  o ^  are  the  observed  fiber  stresses  at  fracture.  When 
the  results  of  Equation  (42)  and  (43),  o'  =  821  psi  and 
m  =  5.94,  were  checked  by  a  Chi- squared  test  the  response  was  at 
the  fifty  percent  level,  indicating  a  very  good  fit. 

For  the  loading  shown  in  Figure  20,  the  beam-column  solution 
is  (ref.  4  ). 
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where  k  =  ~\fp/El.  In  order  to  verify  the  assertion  about  the 
linear  behavior  of  the  beam  up  to  its  cracking  load,  deflections 
were  measured  for  one  beam.  The  cracking  load  was  observed  to  be 
702  pounds  accompanied  by  a  center  deflection  of  .183  inches. 
Substituting  into  Equation  (44)  the  appropriate  values: 


Q  *  702  pounds 
H  =  45  inches 

c  =  7 . 5  inches 
x  =  i/2  *  22.5  inches 


E  =  2.8  x  10^  psi 
I  =  3.26  in4 
P  =  10,000  pounds 
k  =  .0321  in"1 


we  find  y  =  .179  inches,  a  satisfactory  agreement. 


The  bending  stress  in  the  outer  fiber  can  be  found  by  differ¬ 
entiating  Equation  (44) . 


G  =  Ml/21  .  d_  (_  EI  yM) 

(45) 

Edk  sin  kc  ^  r  •  i  i  •  t  /  a 
=  “ZP  ’sin  kl  Q  LSin  kX  +  Sln  k^'x>J 

In  a  manner  analogous  to  that  used  in  developing  Equation  (39)  one 
can  find  the  risk  of  rupture  to  be 
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Using  the  transformation  w  =  x-i/2  and  the  symmetry  about  the 
centerline  this  can  be  written  as 


(47) 


To  normalize  the  integral  we  can  use  z  =  (2/L)w,  leading  to: 


or  in  the  notation  of  Equation  (42), 

a n  m  i  (cr/<x  -l)tn+'*' 

B’7iTIT  </■>  j  — Jh 1 -  dz;  ai0P  (49) 

O  P 

Using  Equation  (45)  to  find  the  o  associated  with  any  value  of 
Q,  Equation  (49)  can  be  evaluated  numerically.  Substitution  into 
Equation  (38)  will  give  the  probability  of  failure  for  any  load  Q. 
The  results  are  shown  in  Figure  23  along  with  the  test  results. 
Bearing  in  mind  that  only  26  prestressed  beam  tests  were  run, 
the  predicted  failure  probabilities  are  fairly  well  borne  out. 
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Figure  23  Probability  of  Failure 


III.  SEGMENTED  COLUMNS 


The  study  of  segmented  columns  in  brittle  materials  is  moti¬ 
vated  in  part  by  the  enormous  potential  of  ceramics  for  resisting 
buckling  and  compressive  fracture  and  in  part  for  the  relationship 
that  exists  between  column  behavior  and  bending  behavior.  There 
are  two  intrinsic  properties  of  segmented  columns  that  demand  spe¬ 
cial  attention.  First,  the  column  has  no  tension  resistance  and, 
consequently,  we  must  modify  the  classical  formulation  of  the  buck¬ 
ling  problem.  Second,  the  imperfect  interface  contract  that  exists 
between  segments  gives  rise  to  statistical  column  behavior.  Fur¬ 
thermore,  rough  interfaces  cause  uneven  loading  across  the  segments 
which  induces  tensile  stresses  in  directions  transverse  to  the  ax¬ 
ial  compression.  The  problems  of  stability,  statistics,  and  strength 
are  studied  in  this  chapter. 


A.  Buckling  -  Perfectly  Flat  Interfaces 


For  the  particular  case  of  a  rectangular  cross  section  with 
perfect  contact  between  segments,  a  static  buckling  analysis  of  a 
segmented  column  will  be  developed.  The  column,  illustrated  in 
Figure  24,  is  assumed  to  be  acted  upon  only  by  the  steady  (nonfol¬ 
lower)  forces  F  a  distance  e  from  the  column's  original  neutral 
axis.  It  is  also  assumed  that  the  number  of  segments  is  very  large 
and  that  the  segment  material  is  entirely  linear  elastic.  Under 
these  assumptions,  solutions  are  found  in  this  section  for  the  sta¬ 
tic  beam  deflection  equations  describing  the  behavior  of  the  seg¬ 
mented  beam-column.  It  is  observed  that  when 


F-Fb 


(50) 


the  deflections  become  unbounded  . 


e 


lim  F 

—  0 


b 


2 

7T 


EIo 


Also  it  is  observed  that 


which  is  the  classic  Euler  buckling  load  for  a  monolithic  column. 
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c)Free  Body  Diagrams  For  Cracked  And  Uncracked  Regions 

Figure  24  Eccentrically  Loaded  Segmented  Column 
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The  development  of  this  solution  will  now  be  presented. 
Referring  to  Figure  24  ,  the  resultant  moment  acting  upon  an 
arbitrary  section  of  the  column  is  found  to  be 

T(x)  =  F  [y(x)  +  e  -  ~y~  f(x)]  (51) 

The  governing  differential  equation  for  the  deflection  curve 
from  elementary  beam  theory  is 

El (x)  =  -  T(x)  (52) 

dxZ 

The  relationships  for  the  crack  penetration  into  a  rectangular 
cross  section  apply  in  this  case,  thus 


f(x)  =0  for  T(x) <  Tc 

f(x)  =  d  -  -6T^  for  T(x)^Tc 


(53) 


where  the  resultant  moment  for  incipient  cracking,  Tc,  is  given  by 


(54) 


The  local  moment  of  inertia  is  given  by 

L  j3 

Kx)  -  IQ  «  for  T(x)  >  Tc 

3 

I(x)  =  18b  [£&!]  for  T(x)  >  T£ 


(55) 


Using  Equation  (51)  and  (54),  Equation  (53)  may  be  expressed  as 


when 

y  < 

T  - e 

then  T(x)  <  Tc 

and  f(x)  =  0 

when 

y  = 

Z  '  e 

then  T(x)  =  T£ 

and  f(x)  =  0 

(56) 

when 

y  > 

*- 

then  T(x)  >  Tc 

and  f(x)  >  0 

Consequently,  the  analysis  will  depend  upon  which  of  the 
following  two  cases  exists:  (1)  entire  column  uncracked,  (2) 
portion  of  column  cracked. 
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Case  (1):  f(x)  =  0 

In  this  case  the  differential  equation,  Equation  (  52),  reduces 
to 

<Ll  +  k2y  =  -  k2e  (57) 

dxz 

where 

k2  =  (58) 

The  boundary  conditions  are  simply 

y  (0)  =  y(L)  =  0  (59) 

The  solution  to  Equation  (57)  and  (59)  is  readily  found  to  be 


y(x)  =  e 


cos  k 


(x-  -j 


cos 


(60) 


Equation 

y  (jr->  - 


(60)  will  be  valid  as  long  as  ymax  <  — g—  -  e.  For  y^ 


~T~ 


F  = 


e,  the  limiting  load  is  found  to  be 


4  E  I 


—  j^cos  ^ 


6e 


1  2 


max 


(61) 


Case  (2) :  f (x)  £  0 

In  this  case  solutions  must  be  found  in  both  the  cracked  and 
uncracked  regions  and  forced  to  match  displacements  and  slopes 
on  the  boundaries  between  the  regions.  Referring  to  Figure  24c, 
x  =  xcL  locates  the  left  boundary  between  the  regions.  Note 
that  T(xcL)  =  Tc  and  y(xQL)  =  -g  -  e.  By  symmetry  xcR  =  L-xcL 
and  thus  i  =  xcR-xcL  =  L-2xcI/ 
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the  solution  found  in 


In  the  uncracked  region,  0  <  x  <  xcL , 
Case  (1)  still  applies,  thus 


y(x)  =  e 


cos  k  (x-  -*-) 


Using  the  fact  that  y(xc^)  =  -  e,  Equation  (62)  yields 


2  -1  ,  d  k  L  N 

E  cos  %I  cos  ~T  > 


In  the  cracked  central  region  &  it  is  readily  shown  that  the 
differential  equation,  Equation  (52),  becomes 


where 


— =  -  h  (y  +  e  -  -|) 

d  x 


TIT 


dv-2 


Equation  (64)  may  be  integrated  once  to  yield 


^Hx^2_2h^y+e“  7  )  1  +  c  (66) 

The  fact  that  it  is  known  a  priori  that  y(xc^)  =  ^  -  e 
enables  c  to  be  determined  at  this  point  by  matching  slopes  and  dis¬ 
placements  in  Equation  (62)  and  (66).  Using  this  value  for  c  and 
Equation  (65)  for  i,  Equation  (66)  becomes 


/dy\  2  _  ,  2  2  t  k  Lv 
(|j)  =  k  e  (cos  -y) 


k  LN  2 


+  2  h  f(y  +  e  -  -j)"1-!- 
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Using  the  fact  that 


'max 


(68) 


and 


Equation  (67)  may  be  solved  for  ymax>  thus 


'max 


-  y  (“2“)  -  T-  -  e  +  d  { 


15 


-  3 


L  d  cos 


3  e 

FT 


-1 


}  (69) 


Inspection  of  Equation  (69)  reveals  that  ymax  becomes 
unbounded  when  the  term  inside  the  braces  approaches  zero,  i.e., 
when 


2 

T 


k  =  cos  (  ££-) 


•1(  6e 
V5d 


or 


4  El 


F  =  Fb  = 


cos-1(^-) 
V5d  J 


(70) 


We  observe  that  in  the  limit  as  the  eccentricity  vanishes, 
approaches  the  Euler  buckling  load  for  monolithic  columns,  thus 


U-  ^ 

e-*-  0 
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B.  Nonflat  Interface  Problem 


1.  Area  and  length  scaling.-  Ideally,  we  would  like  to 
fabricate  segmented  structures  from  segments  with  perfectly  flat 
interfaces;  however,  it  may  not  be  possible  to  obtain  sufficient 
smoothness  on  practical  size  cross  sections.  Even  in  the  case  of 
our  tungsten  carbide  gage  blocks  (one-half  light  band  out  of  flat), 
we  can  still  detect  nonlinearity  in  our  column  response  at  low 
loads.  At  reasonable  temperatures  we  may  be  able  to  approach  mon¬ 
olithic  behavior  through  the  use  of  grout  or  gaskets  between  the 
segments.  At  elevated  temperatures,  however,  these  fillers  may 
have  to  be  dispensed  with.  It  is  for  these  applications  that  we 
will  try  to  improve  our  understanding  of  the  contact  problem. 

From  a  practical  point  of  view,  the  alternative  to  full  scale 
testing  of  segmented  columns  is  to  develop  a  scaling  procedure 
that  will  enable  us  to  predict  the  compression  stress-strain  dia¬ 
grams  for  columns  of  any  length  and  area  from  information  obtained 
from  a  single  segmented  column.  We  begin  our  search  for  such  a 
scaling  law  by  examining  the  effects  of  column  length.  Let  us 
assume  that  all  the  interfaces  have  been  drawn  from  the  same  pop¬ 
ulation  and  that  for  a  specified  axial  stress  there  exists  a  fre¬ 
quency  distribution  for  the  interface  contact  areas.  It  follows, 
then,  that  the  deflection  of  each  segment  will  be  a  random  variable 
which  also  possesses  a  frequency  distribution.  Now,  the  total 
deflection  of  a  multi-segment  column  will  represent  the  sum  of  the 
random  deflections  of  the  constituent  blocks.  Therefore,  the  fre¬ 
quency  distribution  of  the  responses  of  many  nominally  identical 
columns  represents  the  distribution  of  the  sums  of  the  random 
segment  responses.  From  the  Central  Limit  Theorem  of  Statistics 
we  are  assured  that  the  sums  of  random  variables  are  normally  dis¬ 
tributed  regardless  of  the  form  of  the  distribution  for  the  random 
variables  themselves.  Thus,  we  can  hypothesize  that  the  stiffness 
of  segmented  columns  of  a  given  length  and  under  a  specified  load 
are  normally  distributed. 

On  the  basis  that  the  stiffness  of  a  given  size  multisegment 
column  is  normally  distributed,  we  can  proceed  from  the  following 
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theorem  to  describe  the  scaling  rule  for  length.  Theorem  (ref.  5  ): 
"If  x  is  normally  distributed  with  mean  |_i  and  standard  devia¬ 
tion  a  and  a  random  sample  of  size  n  is  drawn,  then  the  sample 

mean  x  will  be  normally  distributed  with  mean  |i  and  standard 

deviation  a//" n."  Evidence  will  be  presented  in  the  next  section 
which  supports  the  hypothesis  of  normality  and  the  scaling  of  the 
scatter. 

It  is  considerably  more  difficult  to  get  a  handle  on  the 
"column  area"  scaling  problem  and,  because  of  a  poorly  executed 
experiment,  our  efforts  in  this  study  can  hardly  shed  any  light 
on  this  matter.  We  do,  however,  have  a  hypothesis  that  we  feel 

is  worth  exploring.  Let  us  examine  the  assumption  that  the  inter¬ 

face  contact  is  controlled  primarily  by  the  highest  asperities  on 
the  surface.  If  the  maximum  asperity  was  measured  on  each  of 
many  nominally  identical  surfaces  we  could  construct  their  fre¬ 
quency  distribution.  The  resulting  frequency  curve  represents 
the  distribution  of  largest  values  in  a  sample  of  size  n  (or 
rather  area  A) .  Methods  for  scaling  such  distributions  to  larger 
areas  are  treated  quite  systematically  by  the  methods  of  extreme 
value  statistics.  Assuming  that  the  stiffness  is  inversely  pro¬ 
portional  to  the  maximum  asperity  heights,  the  stiffness  distri¬ 
bution  F(Et)  might  scale  as  the  distribution  of  smallest  values, 
i.e. , 

Fl(Ec)  -  1  -  [1  -  Fs(Et)]n 

where  the  subscripts  L  and  S  refer  respectively  to  the  large 
and  small  area  columns,  Et  is  the  tangent  modulus,  and 
n  e  A^/Ag.  This  hypothesis  conforms  to  our  past  observations  that 
smaller  area  columns  are  stiffer. 

2 .  Test  results  and  interpretation. - 

a.  Description  of  Experiments:  The  study  of  column  area  and 
length  scaling  requires  a  data  base  of  compressive  stress-strain 
curves  which  we  attempted  to  establish  using  segmented  glass  col- 
lumns .  To  be  sure  that  all  segments  would  be  drawn  from  the  same 
population,  they  were  all  cut  from  the  central  region  of  a  single 
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sheet  of  1/4- inch  window  glass.  Diamond  core  drills  were  used  to 
accurately  produce  fifty  circular  disks  of  each  of  the  diameters, 
1/2,  1,  1-1/2,  2  and  3  inches.  To  fabricate  a  particular  column, 
the  appropriate  size  disks  were  randomly  selected  from  a  rotating 
"lazy  Susan"  on  which  they  were  scattered.  For  ease  in  handling 
and  alignment, the  columns  were  sheathed  in  thin  paper  tubes.  Com- 
pressometers  were  attached  to  tabs  glued  to  the  edges  of  the  top 
and  bottom  disks  in  the  column  as  shown  in  Figure  25.  These  top 
and  bottom  segments  were  reused  for  all  column  tests  of  their  diam¬ 
eter.  Polyethylene  pads  were  placed  on  the  top  and  bottom  of  the 
columns  which  were  then  tested  in  a  Riehle  Universal  Testing 
Machine.  The  load-deflection  diagrams  for  10,  20  and  30  segment 
columns  were  automatically  recorded. 

Seven  typical  load-deflection  diagrams  are  shown  in  Figure  26 
for  two-inch  diameter  columns  with  10  segments.  To  eliminate 
errors  at  low  loads  due  to  backlash  in  the  automatic  plotting 
equipment,  only  unloading  curves  were  considered.  The  highly  indi¬ 
vidual  behavior  associated  with  columns  with  nonflat  interfaces  is 
clearly  illustrated  in  this  figure.  Three  hundred  and  ninety  such 
curves  were  obtained  and  the  hypothesized  statistical  relationship 
which  relates  them  was  investigated  by  selecting  the  particular 
tangent  modulus  associated  with  the  120  psi  compression  level. 

This  stress  was  selected  because  it  fell  close  to  the  knee  in  most 
of  the  curves.  Several  methods  were  studied  for  measuring  the 
slope  of  such  curves  -  two  optical  devices,  a  graphical  technique, 
fitting  a  tangent  by  eye,  and  parabolic  interpolation.  All  of  the 
methods  gave  reasonable  results;  but,  the  latter  method  was  finally 
selected  since  it  produced  the  least  scatter. 

A  simple  computer  program  was  used  to  fit  the  load-deflection 
curves  with  a  parabola,  to  compute  the  slope  at  120  psi,  and  to 
establish  the  mean  and  the  standard  deviation  associated  with  each 
group  of  30  tests  representing  a  given  column  height  and  area. 

The  tendency  to  buckle  precluded  testing  of  the  20  and  30  segment 
columns  for  the  1/2-inch  diameter  disks.  The  cumulative  distribu¬ 
tion  curve  for  the  tangent  modulus  of  each  size  column  was  plotted 
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on  "Normal  Probability  Paper"  and  in  every  case  a  linear  relation¬ 
ship  was  obtained  which  establishes  the  validity  of  our  normality 
hypothesis.  Typical  distribution  curves  are  shown  in  Figure  27 
for  the  two-inch  diameter  column.  The  solid  lines  represent  the 
normal  curves  generated  from  the  computed  mean  and  standard  devia¬ 
tions  of  the  data.  Table  1  tabulates  the  means  and  the  coeffi¬ 
cients  of  variation  (standard  deviation  -5-  mean)  for  each  size 
column. 

b.  Interpretation  of  results:  Referring  all  our  remarks  to 
Table  1  ,  we  first  observe  that  no  consistent  pattern  develops 
vertically  for  either  the  mean  or  the  coefficient  of  variation. 

In  view  of  the  large  scatter,  as  represented  by  the  very  high 
coefficients  of  variation,  this  unexpected  result  may  have  occurred 
by  chance  alone.  On  the  other  hand,  the  data  may  be  in  error. 
Examining  the  parameters  in  the  horizontal  direction  reveals  that 
the  mean  values  consistently  increase  with  increasing  column 
height  and,  except  for  the  three- inch  column,  the  coefficients  of 
variation  uniformly  decrease.  Since  we  cannot  attribute  a  consis¬ 
tent  trend  in  390  observations  to  chance,  we  must  accept  the  pos¬ 
sible  presence  of  a  systematic  error. 

Recalling  that  the  entire  column  and  not  just  the  central 
portion,  was  used  as  the  gage  length,  one  suspects  the  influence 
of  "end  effects."  In  particular,  since  a  properly  run  test  would 
have  a  minimum  of  one  diameter  of  segments  outside  of  the  gage 
length  on  the  top  and  bottom,  it  seems  reasonable  to  suggest  that  the 
Polyethylene  pads  that  were  used  were  too  flexible.  Furthermore, 
this  flexibility  is  of  a  different  character  than  we  would  exper¬ 
ience  with  glass  end  segments.  The  influence  of  such  pads  should 
decrease  with  increasing  column  height  and  this  should  lead  to 
increasing  stiffness.  This  is,  of  course,  exactly  what  we  observe. 
Also,  since  the  same  Polyethylene  pads  and  glass  end  disks  were  used 
for  all  columns  of  equal  area,  we  would  anticipate  that  the  "end 
effect"  would  be  reasonably  constant  for  each  of  the  five  diameters. 
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30  Segments 


Tangent  Modulus  Distribution  Plotted  on  Normal  Probability  Paper 


TABLE  1 


NORMAL  DISTRIBUTION  PARAMETERS  FOR  THE  TANGENT  MODULUS 
OF  SEGMENTED  GLASS  COLUMNS  AT  120  psi  COMPRESSION 


Diameter 

Measure 

Quantity 

Column  Height  (Number  of  Segments) 

10 

20 

30 

1/2" 

mean,  10^  psi 
coef.  var. 

0.92 

28.0  1 

Buckling 

Buckling 

1" 

mean,  10^  psi 
coef.  var. 

1.12 

27.5  7o 

1.33 

9.1  % 

1-1/2" 

mean,  10^  psi 
coef.  var. 

0.63 

21.0  7o 

0.96 

14.1  7, 

1.12 

12.0  % 

2" 

0.70 

28.6  7o 

0.88 

19.1  7, 

1.33 

14.9  7, 

3" 

mean,  10^  psi 
coef.  var. 

0.93 

10.7  7. 

0.93 

10.8  % 

1.12 

13.2  7. 
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Our  contention  that  "end  effects"  influenced  the  behavior  of 
the  segmented  columns  was  briefly  examined  by  testing  an  aluminum 
rod  using  the  identical  procedures  employed  for  the  segmented  glass 
columns.  The  rod  was  10  inches  in  length  with  a  two- inch  diameter 
and  the  initial  portion  of  its  stress-strain  curve  turned  out  to 
be  curvilinear.  The  proper  result  was  recorded  when  the  compresso- 
meter  attachments  were  remote  from  the  ends.  Summarizing  then, 
we  feel  that  our  measurements  of  the  mean  tangent  modulus  do  not 
indicate  the  "pure"  behavior  of  the  segmented  columns  and  that  the 
standard  deviations  are  fair  approximations.  On  this  basis  we 
examined  the  possibility  of  scaling  the  coefficients  of  variations 
of  various  height  columns  and  the  results  are  tabulated  in  Table  2. 
According  to  our  hypothesis,  this  quantity  should  scale  as  1 / \fn. 
Table  2  shows  excellent  agreement  with  this  theory  for  the  1-1/2 
and  2-inch  diameter  columns  and  for  the  average.  This  result  to¬ 
gether  with  the  normality  demonstration  of  Figure  27  tends  to 
support  our  column  height  scaling  hypothesis.  Unfortunately,  no 
conclusions  can  be  drawn  about  the  area  scaling.  Future  experi¬ 
ments  with  segmented  columns  should  incorporate  the  following  im¬ 
provements  : 

(1)  Employ  larger  sample  size  in  view  of  the  large  scatter 
observed . 

(2)  Compressometer  attachments  should  be  kept  in  the  interior 
of  the  column. 

(3)  Digital  output  devices  should  be  utilized. 

C.  Backbone  Column 

It  was  established  in  the  first  phase  of  this  program  that  im¬ 
perfect  contact  between  the  interface  of  a  segmented  column  caused 
transverse  stresses  upon  application  of  axial  loads.  These  trans¬ 
verse  stresses  were  compressive  near  the  interfaces  and  tensile  in 
the  interior  of  the  segment.  Consequently,  it  appeared  that  if  the 
lateral  geometry  of  the  segments  could  be  appropriately  altered,  we 
might  induce  compression  in  the  interior  and  tension  near  the  inter 
faces  of  the  segments.  Thus,  the  two  affects  would  tend  to  cancel 
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TABLE  2 


MEASURED  AND  PREDICTED  VALUES  FOR  THE  COEFFICIENT  OF  VARIATION 
OF  STIFFNESS  FOR  SEGMENTED  GLASS  COLUMNS 


Diameter 

Column  Height  (Number  of  Segments) 

10 

20 

30 

1" 

Measured 

Theory 

27.5% 

27.5% 

9.1% 

19.4% 

8.0% 

15.9% 

1-1/2" 

Measured 

Theory 

21.0% 

21.0% 

14.1% 

14.8% 

12 . 0% 

12 . 1% 

2" 

Measured 

Theory 

28.6% 

28.6% 

19.1% 

.20.2 % 

3" 

Measured 

Theory 

10.7% 

10.7% 

Hi 

Average 

Measured 

Theory 

1 

■ 

■ 

1 
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one  another  and  thereby  give  rise  to  a  segmented  column  with  a 
higher  strength  to  weight  ratio  than  a  prismatic  column.  The 
geometry  selected  for  this  experiment  is  shown  in  Figure  28,  to¬ 
gether  with  a  cylindrical  segment  of  the  same  height  and  interface 
area.  This  latter  specimen  was  used  to  establish  a  group  of  con¬ 
trol  columns . 

Using  the  procedures  outlined  in  Appendix  B,  Hydrostone 
Gypsum  specimens  were  cast  in  each  of  the  shapes  shown  in  Figure  28. 
With  the  simple  setup  shown  in  Figure  29 ,  the  ultimate  compressive 
strength  of  54  three-segment  columns  were  obtained.  The  physical 
and  mechanical  properties  of  these  columns  are  summarized  in 
Table  3 .  It  is  of  considerable  interest  that  a  number  of  columns 
tested  developed  longitudinal  cracks  through  two  or  three  of  the 
segments.  As  shown  in  Figure  29,  the  cracks  did  not  stop  at  the 
interfaces  but  passed  into  the  next  cylinder  as  if  the  column 
were  continuous . 

The  backbone  columns  were  tested  in  the  same  manner  as  the 
cylindrical  columns  as  depicted  in  Figure  30.  As  indicated  in 
Table  3  the  strength-weight  ratio  of  the  backbone  column  is  not 
significantly  higher  than  the  control  column.  Furthermore,  the 
mode  of  failure  casts  doubt  upon  any  possible  superiority  of  the 
backbone  specimen.  Prior  to  ultimate  fracture,  the  lips  or 
flanges  on  the  dogbone  segments  were  stripped  off  leaving  a  pris¬ 
matic  column  of  smaller  diameter.  It  would  then  appear  that  the 
higher  strength-weight  ratio  is  attributable  to  a  size  effect  and 
not  a  geometry  effect. 

When  stress  concentrations  appear  in  a  compressive  field  it 
is  possible  to  achieve  "intelligent  behavior"  from  materials  which 
usually  sustain  no  stress  redistribution  mechanism.  In  addition  to 
the  dogbone  specimen,  another  example  of  such  behavior  was  described 
to  the  author  by  H.  A.  Perry  of  the  Naval  Ordnance  Laboratory. 

Glass  spheres  were  fabricated  from  two  hemispheres  that  were  at¬ 
tached  in  such  a  way  that  a  bead  appeared  around  the  equator  on 
the  inside.  When  the  sphere  was  submerged  in  the  ocean,  the  bead 
was  stripped  off  and  appeared  as  chips  in  the  bottom  of  the  sphere. 
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Figure  28  Example  of  Cylinder  and  Backbone  Specimens 


TABLE  3 

PHYSICAL  PROPERTIES  OF  CYLINDRICAL  AND  BACKBONE  COLUMNS 


Physical  Properties 

Cylindrical 

Backbone 

Interface  Diameter 

3  in 

3  in 

Segment  Height 

3  in 

3  in 

Segment  Weight 

593.6  gm 

333.67  gm 

Number  of  Segments 

3 

3 

Avg.  Ult.  Compressive  Strength 

16989  lbs 

9580  lbs 

Avg. Strength-Weight  Ratio 

Material 

Hydrostone 

Hydrostone 

Number  of  Columns  Tested 

54 

46 

Central  Diameter 

3  in 

3  in 
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As  a  final  observation  concerning  the  backbone  specimen 
we  note  that  this  shape  makes  it  possible  to  apply  a  lateral 
prestress  to  a  column  using  straps  at  a  small  number  of 
locations.  This  idea  is  illustrated  in  Figure  31.  Such  struc¬ 
tures  may  have  advantages  over  jacketed  or  continuously  wound 
columns.  The  bulkier  prestressing  tendons  may  be  less  fragile, 
easier  to  insulate,  or  require  fewer  attachments.  A  similar 
method  of  prestressing  is  briefly  discussed  in  Section  V  for 
the  ogive  shell. 
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IV.  PRESTRESSED  PLATES 


The  similarity  between  prestressed  segmented  beam  behavior 
and  that  of  ductile  bending  made  it  possible  to  successfully  apply 
the  techniques  of  limit  analysis  to  describe  the  ultimate  load 
carrying  capacity  of  prestressed  segmented  beams.  The  extension  of 
these  techniques  to  prestressed  segmented  plates  is  investigated 
in  this  section.  Preliminary  experiments  are  conducted  using  both 
monolithic  and  segmented  circular  plates. 

A.  Segmented  Plates 

If  a  circular  plate  is  subjected  to  a  uniform  radial  pressure 
around  its  periphery,  a  homogeneous  isotropic  state  of  plane  com¬ 
pressive  stress  aQ  is  introduced  into  the  plate.  Any  and  all 
cracks  in  such  a  plate  will  tend  to  close  up,  and  in  particular,  the 
ultimate  bending  resistance  along  such  cracks  will  be  calculated  in 
the  same  manner  used  for  segmented  beams.  Referring  to  Figure  18a, 
the  limiting  moment  per  unit  length  in  a  plate  is  simply 
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(71) 


where  a  is  the  stress  acting  normal  to  the  crack  interface  and 
t  is  the  plate  thickness.  We  shall,  of  course,  take  a  as  the  pre¬ 
stress  Op.  This  limiting  moment  capacity  would .not  be  effected  by 
moments  acting  transverse  to  the  crack  which  suggests  the  applica 
tion  of  the  square  yield  criterion  shown  in  Figure  32  . 
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Figure  32  Square  Yield  Criterion 


76 


We  shall  use  this  criterion  to  analyze  a  circular  plate  of  radius 
R  which  is  simply  supported  on  a  circular  ring  of  radius  r  and  is 
subjected  to  a  central  load  P  brought  onto  the  plate  through  a  blunt 
circular  rod  of  radius  a.  Assuming  the  yield  (or  crack)  pattern 
shown  in  Figure  33,  we  observe  that  the  loading  die  will  ultimately 
contact  the  centerline  of  each  segment  at  only  one  point  a  distance 
of  a  from  the  plate  center.  The  loading  at  such  points  will  be 
P/n  where  n  is  the  number  of  segments.  The  virtual  work  done  by 
these  loads  in  the  assumed  displacement  pattern  is 

n  (  4-  )A  (  ) 


The  energy  dissipated  at  the  yield  lines  is  given  by 


n  R  M  2  sin  -2— 
o  n 


Equating  these  virtual  energies  we  obtain 
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This  load  represents  an  upper  bound  on  the  true  collapse  load  of 
the  plate,  and  consequently,  we  should  choose  from  among  this  class 
of  collapse  mechanisms  the  one  which  gives  the  lowest  load.  This 
occurs  when  n oo,  and  hence. 


p  =  *  %  t2  n_*  “  (73) 

The  true  collapse  load  is  realized  only  when  the  correct  yield  pat¬ 
tern  is  chosen.  In  the  present  case,  symmetry  suggests  that  we 
have  made  the  right  choice. 

We  can  cause  yielding  to  occur  along  a  finite  number  of  radial 
lines  by  strengthening  the  material  between  them.  We  would  expect 
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Deflection  Of  Load  =  A  (-“=*-) 

Rotation  w  *  ^  "  Ae 

n 

Radial  Rotation  Component  =  2ca  sin  -p- 
b  =  r  cos  -~ 


Figure  33  Collapse  Pattern  for  a  Circular  Plate 
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the  capacity  of  the  plate  to  increase  with  such  a  procedure  and 
this  is  exactly  what  Equation  (72)  predicts.  As  we  force  failure 
to  occur  along  fewer  and  fewer  lines,  the  required  strength  of  the 
segments  will  correspondingly  increase.  As  a  rule,  to  avoid  frac¬ 
turing  the  elements  in  a  segmented  component  we  should  select  seg¬ 
ment  geometries  that  approximate  the  true  yield  patterns  for  the 
structure.  An  extensive  treatment  of  limit  analysis  of  plate 
structures  can  be  found  in  Wood  (ref. 6). 

To  establish  the  potential  of  the  proposed  analysis  procedure, 
we  constructed  two  segmented  circular  plates  using  eight  Hydrostone 
plaster  segments  in  each.  The  prestressing  was  accomplished  by 
making  a  double  wrap  of  steel  strapping  about  the  periphery  of  the 
plate  as  shown  in  Figure  34  ,  and  tightening  with  a  standard  band¬ 
ing  tool  until  yielding  occurred  near  the  grip.  Seven  monolithic 
plates  and  two  segmented  plates  were  prestressed  in  this  way  and 
the  resulting  strains  in  six  of  the  plates  were  recorded  by  radial¬ 
ly  positioned  electrical  resistance  foil  strain  gages.  The  strain 
gage  readings  and  measured  loads  for  these  plates  are  tabulated 
in  Table  4  .  The  support  fixture  and  loading  setup  are  shown  in 

Figures  35  and  36  respectively.  The  load  deflection  diagrams  for 
the  segmented  beams  are  shown  in  Figure  37  where  we  observe  well 
defined  horizontal  regions.  Upon  unloading,  we  obtained  complete 
deflection  recovery  with  only  occasional  chipping  at  the  segment 
edges . 

As  evidenced  from  Table  4  ,  the  straps  did  not  apply  a  uniform 
radial  prestress;  however,  since  care  was  taken  to  tighten  the  straps 
in  the  same  way  for  all  cases,  it  is  felt  that  the  average  prestrains 
in  the  various  plates  were  about  the  same.  On  this  basis,  the  pre¬ 
strain  was  taken  as  the  average  at  sixteen  gage  readings,  i.e., 
49.4x10  ^  inch/inch.  Using  this  value  together  with  the  plate  pro¬ 
perties  tabulated  in  Table  5  ,  Equation  (72)  predicts  a  yield  load 
of  P  =  309  lb.  This  value  differs  from  the  measured  values  of 
270  lb  and  288  lb  by  14.4  percent  and  7.3  percent  respectively. 
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TABLE  4 


PRESTRESSED  MONOLITHIC  AND  SEGMENTED  CIRCULAR  PLATE  STRENGTHS 


Plate  Number 

Strain 

x  10"  6 

in . / in . 

Yield  Load 

Ultimate  Load 

a 

b 

c 

(lb) 

(lb) 

Monolithic  1 

defec  t. 

50 

45 

400 

540 

Monolithic  2 

70 

55 

65 

345 

582 

Monolithic  3 

30 

50 

55 

370 

520 

Monolithic  4 

40 

45 

55 

410 

624 

Monolithic  5 

60 

70 

defect . 

330 

706 

Monolithic  6 

n 

d  gages 

300 

535 

Monolithic  7 

n< 

d  gages 

350 

430 

Segmented'  1 

50 

40 

270 

— 

Segmented  2 

n 

_ 

3  gages 

_ 

288 

— 

—  6 

Average  Strain  Gage  Reading:  49.4x10"  in. /in. 
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Figure  35  Ring  Support  for  Circular  Plates 
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Deflection 

Figure  37  Central  Load  -  Central  Deflection  Diagrams 
for  Prestressed  Segmented  Circular  Plates 


TABLE  5 


PHYSICAL  PROPERTIES  OF  HYDROSTONE  PLATES 


Plate  radius 

R  = 

7.5  in. 

Support  ring  radius 

r  = 

6.75  in. 

Central  load  die  radius 

a  = 

0.906  in. 

Plate  thickness 

t  = 

5/8  in. 

Average  plate  strain 

£p  * 

49.4x10  ^  in. /in. 

Modulus  of  elasticity  of  Hydrostone 

E  = 

2.79x10^  psi 

Poisson's  Ratio  for  Hydrostone 

v  = 

1/4 

Average  prestress  level 

a  = 
P 

Epl?v=  184  Psi 

Number  of  segmented  beam  elements 

=  8 
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B.  Monolithic  Plates 


When  a  prestressed  monolithic  plate  is  slowly  loaded,  conven¬ 
tional  elastic  behavior  is  experienced  until  the  net  tensile  stress 
at  some  point  exceeds  the  material  strength  and  a  crack  develops. 

In  a  constant  strain  rate  machine,  the  load  would  fall  off  abruptly 
and  then  increase  again  as  the  strain  continued  to  increase.  This 
behavior  is  depicted  in  the  central  load-central  deflection  diagram 
shown  in  Figure  38  for  the  second  monolithic  plate  referenced  in 
Table  4.  As  we  see,  other  cracks  continue  to  form  until  the  strength 
of  the  surviving  material  is  sufficient  to  force  unconstrained  yield¬ 
ing  in  the  existing  crack  pattern.  The  radial  crack  pattern  (artifi¬ 
cially  darkened)  shown  in  the  top  plate  in  Figure  34  is  exactly  what 
our  limit  analysis  theory  anticipates.  The  visible  crushing  at  the 
center  of  this  early  plate  test  was  caused  by  a  steel  ball  that  was 
used  originally  to  load  the  plate.  Subsequent  tests  employed  a  cir¬ 
cular  die  to  distribute  the  central  load  over  a  greater  area. 

Seven  central  load  tests  were  conducted  with  monolithic  plates 
and  their  yield  and  ultimate  strengths  are  recorded  in  Table  4  . 

The  average  ultimate  strength  for  these  members  is  562  pounds  which 
represents  a  considerable  increase  in  the  average  strength  of  mono¬ 
lithic  plates  without  prestressing,  328  pounds.  The  distribution 
curve  for  the  strength  of  these  unprestressed  control  plates  is 
shown  in  Figure  39  where' we  observe  a  considerable  spread  in  the 
data.  This  implies  that  very  low  strength  values  will  be  present 
in  a  population  of  even  a  few  hundred.  Consequently,  very  "low 
strength"  operating  levels  must  be  used  to  obtain  reasonable  relia¬ 
bility.  For  the  prestressed  plates,  on  the  other  hand,  there  is  a 
built  in  fail-safe  mechanism.  Although  the  ultimate  load  is  statis¬ 
tical  and  may  be  subject  to  wide  variability,  the  yield  load  is 
bounded  from  below.  We  observe  that  the  weakest  plate  is  achieved 
when  an  infinite  number  of  radial  cracks  develop  in  which  case 
their  yield  load  is  computed  from  Equation  (73) .  Real  plates  will 
crack  in  a  finite  number  of  places  and  will  therefore  be  stronger. 

For  our  prestressed  monolithic  plates,  Equation  73  predicts  a  yield 
load  of  289  pounds  and  we  observe  from  Table  4  that  all  of  the 
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(lbs) 


inches 

Figure  38  Central  Load  -  Central  Deflection  Diagram  of  a  Prestressed 
Monolithic  Circular  Hydrostone  Plaster  Plate 


yield  strengths  reported  are  higher  than  this  value;  the  average 
load  is  358  pounds . 

We  should  point  out  that  our  two  mathematical  models  for 
beams  assume  an  infinite  number  of  segments  and  should  therefore 
predict  a  lower  bound  on  the  behavior  of  monolithic  prestressed 
beams.  In  the  exceptional  case  of  a  beam-column;  however,  we  can¬ 
not  make  this  statement  since  the  resistance  at  a  section  can  be 
lowered  by  virtue  of  the  beam  deflection.  The  distance  from  the 
outer  compressive  beam  fibers  to  the  centroid  of  the  tendons  de¬ 
creases  as  the  beam  deflection  increases.  As  a  matter  of  fact, 
a  point  of  instability  is  finally  reached  and  the  prestressing 
force  participates  in  the  catastrophic  destruction  of  the  beam. 
Finally,  since  a  monolithic  segment  sustains  no  separation,  it  is 
stiffer  than  an  equivalent  length  of  several  segments.  Consequently, 
beams  constructed  with  a  finite  number  of  elements  will  be  stiffer 
than  predicted  from  an  infinite  number  of  elements. 
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V.  PRELIMINARY  INVESTIGATIONS  OF  PRESTRESSED  SHELLS 

Fundamental  studies  of  prestressed  elements  demand  as  a  pre¬ 
requisite  that  methods  be  available  for  applying,  maintaining,  and 
monitoring  prestressing  forces.  In  the  case  of  shell  structures, 
this  may  require  a  considerable  amount  of  innovation  and,  perhaps, 
the  elaborate  development  of  special  techniques.  Since  such 
efforts  are  beyond  the  scope  of  our  present  endeavors,  a  brief 
study  has  been  undertaken  to  examine  two  conventional  approaches 
to  the  problem  of  prestressing  an  ogive  shell  and  a  cylindrical 
shell. 


A.  Cylindrical  Shell 


For  problems  which  require  the  continuous  application  of  pre¬ 
stressing  force  over  large  areas,  it  is  sometimes  possible  to  utilize 
the  technique  of  overwinding.  Here,  we  visualize  using  a  brittle 
component  in  place  of  a  mandrel.  Continuous  filaments  under  high 
tension  are  then  wound  onto  the  components,  and  in  this  way,  a  jacket 
is  formed  permanently  over  the  monolithic  or  segmented  element  which 
may  prevent  the  leakage  of  liquids  or  gases  and  which  could  provide 
an  energy  absorbing  layer  that  would  protect  the  component  from 
local  impacts . 

To  investigate  the  potential  of  this  technique,  a  three  inch 
diameter  alumina  cylinder  with  a  0.2  inch  wall  thickness  was  over¬ 
wrapped  with  preimpregnated  20  end  glass  roving.  An  electrical 
resistance  strain  gage  was  mounted  on  the  inner  wall  of  the  cylin¬ 
der  to  monitor  the  induced  stress  caused  by  the  overwinding.  Using 
an  overwinding  tension  of  0.325  lbs/end  to  produce  a  tension  force 
of  65  lb/in.  along  the  cylinder  axis,  we  obtained  the  relationship 
shown  in  Figure  40  between  induced  prestress  and  fiberglass  layers. 


The  relationship  between  the  induced  prestress  in  a  cylinder 
and  the  number  of  layers  of  overwound  filaments  was  obtained 
where  account  is  taken  of  the  changes  in  filament  stress  that  occur 


when  additional  filament  layers  are  applied, 
is  given  by 


The  hoop  stress  a 
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Theory 


of  Overwound  Fiberglass 


where  5  is  the  layer  thickness,  t  is  the  wall  thickness  of 
the  cylinder,  r  is  the  outside  radius  of  the  cylinder,  T  is  the 
overwind  tension  per  unit  length  in  the  direction  of  the  axis  of 
the  cylinder,  and  E-^  and  E2  are  the  modulus  of  elasticity  of 
the  cylinder  material  and  the  filament  material  respectively. 

This  theory  assumes  that  we  have  perfect  packing  of  rectangular 
fibers;  however,  for  circular  fibers  without  nesting  we  can  modi¬ 
fy  the  theory  by  replacing  E2  by  7TE2M.  Using  the  following, 
data,  the  theory  is  plotted  in  Figure  : 

T  =  65  lb/in.  t  =  3/32  in. 

5  =  0.005  in.  E^  =  40x10^  psi 

r  =  1.5  in.  E2  =  8x10^  psi. 

Considering  the  preliminary  nature  of  the  experimental  setup  used, 
the  agreement  between  the  measured  and  predicted  values  is  quite 
satisfactory.  On  this  basis,  calculations  performed  for  ultra 
high  strength  overwinding  materials  such  as  high  carbon  steel  wire 
and  "S"  fiberglass  single  end  indicate  that  induced  stresses  ex¬ 
ceeding  1000  psi/layer  for  a  one  inch  wall  thickness  specimen  are 
not  unreasonable.  This  implies  that  for  practical  purposes  any 
desired  level  of  prestress  can  be  obtained  provided  that  winding 
on  a  geodesic  path  is  possible. 

B .  Ogive  Shell 

If  an  axial  compressive  load  is  applied  to  the  apex  of  an  ogive 
shell,  the  membrane  stresses  in  the  shell  wall  are  compressive  in 
every  direction.  We  see  from  this  example  that  an  effective  state 
of  initial  stress  can  be  accomplished  by  applying  a  prestressing 
force  at  a  discrete  point.  Because  the  application  and  monitoring 
of  such  a  force  is  exceptionally  simple,  the  fabrication  of  ogive 
shells  was  attempted  during  this  phase  of  the  program. 

Figure  41  illustrates  one  of  several  segmented  Hydrostone 
plaster  shells  which  were  produced  by  cutting  up  monolithic  shells 
before  they  completely  cured.  The  interfaces  of  these  shells  tend 
to  close  up  under  compression;  however,  the  thin  sections  near  the 
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Figure  41  Segmented  Ogive  Shell  Prestressed  at  it's  Apex 


apex  are  breaking  under  load.  We  attribute  these  failures  to 
extremely  rough  interface  conditions.  Further  studies  with  these 
shells  should  address  themselves  in  this  problem  in  addition  to 
investigations  of  prestressed  monolithic  behavior. 
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APPENDIX  A 


I-BEAM  COMPUTER  PROGRAM 

In  this  appendix,  we  shall  present  the  computer  program  which 
was  developed  to  perform  the  analysis  of  prestressed  segmented 
I-beams  with  multiple  elastic  tendons.  The  program  consists  of  a 
main  program  plus  ten  function  subprograms.  The  main  program  es¬ 
sentially  consists  of  those  equations  presented  in  the  section  on 
general  relationships.  Those  relationships  which  depend  explicitly 
upon  the  particular  cross  section  being  investigated  have  been  iso¬ 
lated  into  individual  function  subprograms..  The  applied  bending 
moment  distribution,  in  normalized  form,  has  also  been  isolated 
in  its  own  function  subprogram.  Hence,  only  the  appropriate  sub¬ 
program^)  need  be  changed  to  permit  the  investigation  of  different 
cross  sections  and/or  applied  loadings. 

Immediately  following  this  introduction  is  a  glossary  of  some 
of  the  important  variables  as  they  appear  in  the  program  along  with 
their  relationship  to  the  variables  used  in  the  presentation  of  the 
analysis  elsewhere  in  this  report.  Next  after  the  glossary  is  a 
listing  of  the  computer  program  as  written  in  Fortran  II  for  the 
IBM  7094.  Finally,  a  sample  input  and  corresponding  output  are 
presented. 

The  output  has  been  designed  to  present  as  much  relevant  infor¬ 
mation  as  possible  in  a  meaningful  manner.  The  first  page  of  out¬ 
put  contains  all  the  input  data  plus  the  results  of  a  few  preliminary 
calculations.  Next  is  presented  the  table  of  W  vs.  f.  The  first 
differences  in  W  are  also  given  to  illustrate  whether  the  function 
is  approximately  linear  between  entries  in  the  table  as  assumed  by 
the  linear  interpolation  scheme  in  the  function  subprogram  FCP(W) . 

The  remainder  of  the  output  consists  of  evaluating  deflections, 
end  rotations,  etc.  at  the  various  selected  values  of  the  applied 
loading.  The  first  value  of  the  applied  load  is  always  zero  and  is 
followed  by  an  intermediate  value  which  leads  to  the  value  corres¬ 
ponding  to  incipient  cracking-.  The  program  then  truncates  the 
cracking  load  to  a  specified  number  of  figures  and  proceeds  by  adding 
specified  increments  in  applied  load  until  either  the  load  exceeds 
the  specified  maximum  or  one  of  the  counters  exceeds  its  specified 
maximum.  The  final  item  in  the  output  is  a  table  of  applied  load 
versus  central  deflection,  left  end  rotation,  and  right  end  rotation. 
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Partial  listing  of  variables  appearing  in  computer  program: 


AL  .  applied  load,  e.g.  force,  moment,  etc. 

ALDMBM .  applied  load  divided  by  maximum  bending  moment, 

eg.  P/ (PL,/ 4) 

ALMAX . value  of  applied  load  which  is  not  to  be  exceeded 

AT  .  area  of  each  tendon,  At 

CP  .  crack  penetration, f 

CPXD  . crack  penetration  f  at  x^ 

DAL  .  increment  in  applied  load 

DCENT  .  distance  q  from  bottom  of  section  to  center  of 

gravity 

DELTA .  deflection  A  at  xA 

DIM1  =  B . width  of  flange  for  I-beam,  b 

DIM2  =  TW  .  thickness  of  web  for  I-beam,  t 

w 

DIM3  =  D  .  depth  of  web  for  I-beam,  d 

DIM4  =  TF  .  thickness  of  flange  for  I-beam,  t^ 

DIM5,  etc.  =  0..  not  used  for  I-beams 

DISC(I)  ........  distance  from  bottom  of  section  to  i  discontinuity 

DTOT  .  overall  depth,  dtota]_,  of  cross  section 

EB  .  elastic  modulus  of  beam,  E 

EC  .  effective  initial  eccentricity  of  preload,  e 

■f-Ti 

ECT(I)  .  eccentricity  of  i  tendon,  e^ 

ET . . .  elastic  modulus  of  tendon,  Et 

F  . .  total  preload  exerted  by  tendons,  F 


F0 . . .  initial  preload  exerted  by  tendons,  FQ 

FT(I)  .  force  in  the  it!l  tendon,  F^ 

FTO(I)  . .  initial  force  in  it'ri  tendon,  F0^ 
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IK .  number  of  equally  spaced  points  on  beam  at  which 

deflections  are  to  be  found 

IKK .  if  IKK  =  2,  then  applied  loading  is  symmetrical 

about  x  =  L/2 

K  .  number  of  increments  used  in  numerical  integrations 

KFUDMl . maximum  number  of  applied  load  increment  halvings 

KFUDM2  .  maximum  number  of  occurrences  of  questionable  output 

KSMAX .  maximum  number  of  steps /loop  in  F  and  FL, 

directions  in  the  solution  of  tendon  equations 

KTRYM . maximum  number  of  loops  in  the  solution  of  tendon 

equations 

NIR(I)  . number  of  entries  in  table  of  W  vs.  f  between 

f  =  CP  =  DISC (1-1)  and  DISC (I) 

NOREG  .  total  number  of  regions  cross  section  is  divided 

into  depthwise  for  table  of  W  vs.  f  (three  for 
I-beam) 

NSF  .  number  of  significant  figures  for  truncation  of 

cracking  load 

NT  .  total  number  of  tendons,  n 

PERMAX . maximum  allowable  percent  in- out  errot  in  the 

solution  of  the  tendon  equations 

T  .  resultant  bending  moment  applied  to  beam,  T 

TABAL(I)  .  i  applied  load  for  table  at  end  of  output 

TABCD(I)  .  it41  central  deflection  for  table  at  end  of  output 

4-t-. 

TABCP(I)  .  i  crack  penetration  for  table  of  W  vs.  f 

i-|_ 

TABTL(I)  .  1  left  end  rotation  for  table  at  end  of  output 

j_l_ 

TABTR(I)  .  iL  right  end  rotation  for  table  at  end  of  output 

TABW(I)  .  it41  value  of  W  for  table  of  W  vs.  f 

THETAL  .  left  end  rotation,  9^ 

THETAO  .  initial  end  rotation,  0q 

THETAR .  right  end  rotation,  0_ 

_  ft 

V  . maximum  value  of  W  along  the  beam 
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w 


variable,  W  =  (M-tO/F,  appearing  in  crack 
penetration  equation 

WCRACK  .  value  of  W  corresponding  to  incipient  cracking 

X  .  coordinate,  x,  measured  along  the  length  of  the 

beam 

XD  .  location,  x^,  where  deflection  is  being  determined 

YL  . .  lenth  of  segmented  beam,  L 

YLT  .  length  of  tendons,  Lt 

YMB  . maximum  applied  bending  moment, 

YMT  .  tendon  stiffness  bending  moment, 
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-n  -n  ti  -n 


PROGRAM  LISTING 


CMAIN  PROGRAM  -  PRESTRESSED  SEGMENTED  I  BEAM  WITH  ELASTIC  TENDONS,  PCH 
DIMENSION  HOLLER! 12)»DIFX(200) 

DIMENSION  FTO( 10),ECT( 10) , FT! 10) , S I GMA ( 10) , TABAL ( 100 ) , TABTL ( 100 ) , 
IT  ABTR{ 100), DISC! 10) ,NIR( 10 ) , TABCP < 200 ) , TABW! 200 ) , TABCD ( 100 ) 

COMMON  YMB,YMT,F,FO, YL,EB,XD,WCRACK,NI , IC , IBMD , I XS , DTOT ,DCENT,D I  Ml 
1,DIM2,DIM3,DIM4,DIM5,DIM6,DIM7,DIM8,DIM9,DIM10,DI SC, TABCP, TABW, EC 
FDI 
FTLI 
FTR  I 
FBC I 

45  WRITE  OUTPUT  TAPE  6,50 

50  FORMAT! 55H IPRESTRESSED  SEGMENTED  I  BEAM  WITH  ELASTIC  TENDONS,  PCH) 
WRITE  OUTPUT  TAPE  6,55 

55  FORMAT! 119H0DIMENSI0NS  -  LENGTHS  ARE  IN  INCHES,  LOADS  ARE  IN  KIPS, 

I  BENDING  MOMENTS  ARE  IN  INCH-KIPS,  STRESSES  +  MODULI  ARE  IN  KSI  ) 

100  READ  INPUT  TAPE  5, 101, HOLLER 

101  FORMAT! 12A6) 

READ  INPUT  TAPE  5, 105 , YL , E B , YLT, ET , AT , PERMAX, K , NT , NOREG 
105  F0RMAT!6F10. 0,315) 

READ  INPUT  TAPE  5,110, tFTO! I),ECT!  I), 1  =  1, NT) 

110  FORMAT ( 2F10.0) 

READ  INPUT  TAPE  5, 1 15 , DTOT, DCENT , D I  Ml ,D IM2 ,D I  M3 , DI M4 
READ  INPUT  TAPE  5 , 1 1 5, D IM 5, D IM6 , DI M7 , D I M8 , DI M9 ,D I M10 
115  FORMAT ( 6  F 1 0 . 0 ) 

READ  INPUT  TAPE  5,120, ! 0 1  SC ! I) ,NIR! I ) ,1=1 .NOREG ) 

120  FORMAT  IFIO.6,15) 

READ  INPUT  TAPE  5 , 1 27 , I K , I KK , NSF , ALDMBM, ALMAX , DAL 

127  FORMAT! 315, 3F10.0) 

READ  INPUT  TAPE  5, 128, KTRYM ,KSMAX» KFUDM1 .KFUDM2 

128  FORMAT (415) 

WRITE  OUTPUT  TAPE  6, 130, HOLLER 
130  FORMAT! 1H0,12A6) 

WRITE  OUTPUT  TAPE  6, 1 35 , YL , EB , YL T, E T, AT 
135  FORMAT! 15H0BEAM  LENGTH  =  ,F8.4,12H  8EAM  MOD  =  ,F10.2,17H  TENDON  LE 
1NGTH  =  , F8 .4, 14H  TENDON  MOD  =  ,F10.2,19H  AREA  PER  TENDON  =  ,F7.5> 
WRITE  OUTPUT  TAPE  6, 140, PERMAX, K 

140  FORMAT! 75H0PERCENT  IN-OUT  DIFFERENCE  IN  PRESTRESS  AND  TENDON  MOMEN 
IT  IS  LESS  THAN  OR  =,F6.4,33H  NO.  INCREMENTS  IN  INTEGRATIONS  =,I4) 
WRITE  OUTPUT  TAPE  6 , 14 5 , NT , NOR EG 
145  FORMAT! 18H0NQ.  OF  TENDONS  =  ,I2,18H  NO.  OF  REGIONS  =  ,12) 

WRITE  OUTPUT  TAPE  6 , 150, DTOT , DCENT , D I M 1 , DI M2 ,D I M3 , DI M4 
150  FORMAT ( 15HOTOTAL  DEPTH  =  ,F7.4,20H  DIST  TO  CENTROID  =  ,F7.4,8H  DIM 

II  =  ,F7.4,8H  DIM2  =  ,F7.4,8H  DIM3  =  ,F7.4,8H  DIM4  =  ,F7.4) 

WRITE  OUTPUT  TAPE  6 , 1 55, D I M5 , D I M 6, D I M7 , DI M8 , DI M9 ,D I  Ml 0 

155  FORMAT ( 8H0DIM5  *  ,F7.4,8H  DIM6  =  ,F7.4,8H  DIM7  =  ,F7.4,8H  DIM8  =  , 
1F7.4.8H  DIM9  =  ,F7.4,9H  DIM10  =  ,F7.4) 

WRITE  OUTPUT  TAPE  6 , 1 60, ( I , D I  SC ( I ) , I ,N IR ( I ) , I = 1 , NOREG) 

160  FORMAT ! 6H  DISC!, 12, 5H  )  =  ,F7.4,5H  NIR(,I2,5H  )  =  ,13) 

SUMECT  =  0. 
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DO  165  1=1, NT 
165  SUMECT=SUMECT+ECT( I ) 

ERRM AX =PERM AX/IOO. 

IXS  =  0 
I 8MD=0 

WRITE  =  FM(O.O) 

WRITE  =  FIIO.OI 

IXS  =  1 

I BMD=1 

XNT=NT 

R=0. 

RM=0. 

DO  170  1=1, NT 
FT ( I )  =  FTO(  I) 

R=R*FTO(I) 

170  RM=RM«-FTO(  I)*ECT(I) 

FOR 

EORM/R 

F=F0 

WRITE  OUTPUT  TAPE  6 , 17 5 , ALDMBM , ALMA X , DAL 
175  F0RMAT(35H0( APPLIED  L0AD)/(MAX  BEND  MOMENT)  =,F10.6,20H  MAX  APPLIE 
ID  LOAD  =  , F10.6, 29H  INCREMENT  IN  APPLIED  LOAD  =  ,Fi0.6) 

WRITE  OUTPUT  TAPE  6, 180, F, EC 

180  FORMAT (27H0INITIAL  PRELOAD  IN  8EAM  =  ,F10.6,35H  ECCENTRICITY  OF  IN 
1ITIAL  PRELOAD  =  ,F10.6) 

WRITE  OUTPUT  TAPE  6,185, IK, NSF 

185  FORMAT ( 18H0 DEFLECT  IONS  AT  1/,I2,16H  POINTS  ON  BE AM , , 28H  ROUND  OFF 
1CRACKI NG  LOAD  TO  ,I1,20H  SIGNIFICANT  FIGURES  ) 

I F ( IKK-2 ) 190,187, 190 

187  WRITE  OUTPUT  TAPE  6,188 

188  FORM AT (40H0AP PLIED  BENDING  MOM  DIST  IS  SYMMETRICAL  ) 

190  WRITE  OUTPUT  TAPE  6 , 195 , K TRYM , KSMA X 

195  FORMAT (70H0F0R  DETERMINATION  OF  PRELOAD  AND  TENDON  MOMENT  -  MAX  NO 
1.  STEPS/LOOP  =  , 12, 16H  MAX  NO.  LOOPS  =  ,I2,21H  MAX  NO.  APPLIED  LOA 
2D  ) 

WRITE  OUTPUT  TAPE  6 , 196, KFUDMl , KFUDM2 

196  FORMAT ( 21H0 INCREMENT  HALVINGS  =,I2,44H  MAX  NO.  OCCURENCES  OF  QUEST 
1I0NABLE  OUTPUT  =,I2) 

C  TABLE  OF  TABW  VS  TABCP  (200  THRU  299) 

200  N I R1P1 =N I R ( 1 ) + 1 

DO  210  K2= 1, N IR IP  1 

TABCP( K2)=FL0ATF(K2-i)*DISC( 1 ) /FLO A TF ( N I R ( 1 ) ) 

X1=TABCP(K2) 

210  TABW(K2)  =  EC+X1-DCENT+FXB(X1)+FI(XI  )/(FA(Xl >*FXB(Xl>  > 

WCRACK=TABW( 1) 

N I RT=N IR( 1 ) 

I F ( NOREG-2  )  240,215,215 
215  DO  230  K1=2,N0REG 
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X2  =  DISC(K1 )— 0ISC(K1— 1) 

NNN=NI R( K 1 ) +1 
DO  220  K2= 1 , NNN 
K3  =  NIRT  +K2 

TABCP{K3)=DISC<KI-1)+FL0ATF!K2-1)*X2/FL0ATFINIR<K1) ) 
X1=TABCP!K3) 

T  ABW( K3 )® EC+X1— DCENT+FXB ( XI ) +F I C  XI ) /C  FA  I XI  )>FXB( XI I ) 

220  CONTINUE 

NIRT=NIRT+NIR(K 1) 

230  CONTINUE 
240  N I  =  N IRT 

WMAX=TABW(NI+1) 

245  WRITE  OUTPUT  TAPE  6,250, WCRACK 
250  FORMAT! 9H1WCRACK  =,F10.6) 

WRITE  OUTPUT  TAPE  6,255 

255  FORMAT ( 1H0, 13X.6H  TABCP,14X,5H  TABW,12X,11H  FIRST  OIFF  ) 
DI FX ( L ) =0.0 
N I  P2=N I +2 
DO  270  K4=2,NIP2 
DIFX(K4)=TABW(K4)-TABW!K4-I» 

K  5  =  K  4—  1 

WRITE  OUTPUT  TAPE  6 , 260,  K5»  TABCP  I  K5  )-,  TABW1 K5 ) *DIFX  !  K5 ) 
260  FORMAT! I5.3F20.8) 

270  CONTINUE 

DO  280  1=1, NT 
280  FT (  I  )=  FTO!  I) 

CTEN=YLT/I AT*ET) 

V=0. 

YMT=0. 

MM=0 

YMB=0. 

X  I K= IK 
I JK=  IK/ IKK 
KK=0 
KFUD=0 
X=0. 

TABTL! 1)=FINT(FTLI,K) 

TABTR! 1 )  =  F I NT (FTRI»K) 

TABALI l)=0. 

THETAO=TABTL( 1) 

THET  Al_  =  THET  AO 
THETAR=THETAO 
WRITE  OUTPUT  TAPE  6,290 
290  FORMAT (1H1) 

WRITE  OUTPUT  TAPE  6,295, THETAO 
295  FORMAT! 10H  THETAO  =  ,F12.8) 

GO  TO  400 

C  ROUTINE  FOR  SELECTING  APPLIED  LOAD  ! 300  THRU  399) 

300  KK=KK+1 
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c 


c 


310 

320 

330 

340 

345 

350 

355 

400 

410 


411 

412 

413 

414 


415 

418 


419 

425 


I F ( KK-2  )  310,320,330 
YMB=WCRACK*F0 
YMB0LD=YM8 
GO  TO  400 

YMB=( ( FO»WCRACK ) **2  )/( ( 2.*F0-F ) *WCRACK-YMT ) 

ID=XFI XF ( . 43429448*L0GF { YMB* ALDMBM ) )-NSF+l 
ALO=(lO»**IO)*FLOATF(XFIXF( YMB* ALDM8M/ ( 10. ** I D ) ) ) 


SFYMB= ALD/ ALOMBM 
GO  TO  400 

YMB=SFYMB+DAL#FLOATF< KK-2) /ALDMBM 


AL= ALD+DAL»FLOAT  F { KK-2 ) 

IF ( AL- ALM AX >400,400,340 
WRITE  OUTPUT  TAPE  6,345 

FORMAT ( 15H1  APPLIED  L0AD,8X,12H  CENTER  DEFL, iOX, 11H  THETA  LEFT ,8X 
1, 12H  THETA  RIGHT) 

WRITE  OUTPUT  TAPE  6,350, ( TABAL ( N) , TABCD l N) ,TABTL l N> , TABTRI N) ,N= 1 ,M 
1M) 

FORMAT ( 4F20 • 8 ) 

WRITE  OUTPUT  TAPE  6,355 
FORMAT  ( 1H 1 ) 

GO  TO  45 

ROUTINE  FOR  DETERMINING  PRESTRESS  FORCE  AND  MOMENT  (400  THRU  499) 
AL-YMB*ALDMBM 


THETAL=FINT(FTLI»K) 

THET  AR  =  FINT ( FTR I , K ) 

IF(ET*FLOATF(KK) )  410,500,410 

ET=0.  FOR  ZERO  STIFFNESS  TENDONS 

CYMT-Y  MTOLD 

YMT=YMT  OLD*YMB/YMBOLD 

F=FOLD* ( YMB-YMT ) /( YM BOLD- YM TOLD ) 


DELF=(F-F0LD)/2. 

DELYMT  = ( YMT-YMTQLD ) /2. 


KKFU0=0 

IF ( ABSF(DELF/F)-ERRMAX >41 1,411,412 
0ELF=2.*ERRMAX*F 

IF (ABSF( DELYMT )-ABSF I YMT*ERRMA X ) >413,413,414 

0ELYMT-2.*ERRMAX*YMT 

SF  =  DEL  F 

S YMT=DELYMT 

KS=0 

WRITE  OUTPUT  TAPE  6,415,AL 
FORMAT ( 15H0APPLIED  LOAD  =,F12.8) 


KTRY*0 

KALT=0 

WRITE  OUTPUT  TAPE  6,419,YMT 
FORMAT! 14H  HOLDING  YMT  =,F12.8) 


DFOLD=F 

FIN=F 

IF(KTRY-KTRYM) 426, 426,461 
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426  WL=f(YMB»FM(0-  )-YMTI/F 
CPL=FCP(  WL  ) 

WR=(YMB*FM(YL)-YMT)/F 
CPR=FCP( WR ) 

THETAL=FINT(FTLI,K) 

I F I IKK-2 )  428,427,428 

C  I KK=2  FOR  SYMMETRICAL  LOADING,  FM ( X ) =FM( L-X) 

427  THET AR=THETAL 
GO  TO  429 

428  THET AR=  FINT(FTRI,K) 

429  BCOMP=FINT(FBCI»K) 

YMT  OUT  =  0. 

Z1=(SUMECT+XNT*(CPL— DCENT+FXBtCPL) > )*< THETAL-THETAO) 

Z  2= ( SUMECT  +XNT  *(CPR-DCENT+FXB(CPR) )) * ( THETAR-THETAO) 
F0UT=F0+(Z1+Z2)/(CTEN*XNT*BC0MP  ) 

DO  430  1=1, NT 

ZZ 1=  (  ECT  (  I  )+CPL-DCENT*-FXB(CPL)  >*<  THETAL-THETAO) 

ZZ  2=  (  ECT(  I  )<-CPR-DCENT  +  FXBICPR)  )  *  I  THE  TAR- THETA 0 ) 

FT ( I )=  FTO ( I)+(ZZ1*ZZ2-(F-F0)*BC0MP)/CTEN 
4  30  YM TOUT  =  YMT CUT *FT  ( I) * ( ECT ( I > -EC ) 

DF=FOUT-F IN 
435  VIN=(YMB-YMT)/FIN 

VOUT=( YM8— YMTOUTJ/FOUT 

WRITE  OUTPUT  TAPE  6 , 438 , F IN , V I N , FOUT , DF , YM TOUT , VOUT 
438  FORMAT ( 6H  FIN  =,F12.8,6H  VIN  =,F10.8,7H  FOUT  =,F12.8,5H  DF  =,F12.8 
1 , 9H  YMTOUT  =,F12.8,7H  VOUT  =,F10.8) 

IF( ABSF( (FOUT-FIN)/FIN)-ERRMAX)  450,450,440 
440  IF(DFOLD*DF)445,450,442 

442  I F  < ABSF(DF0LD)-A8SF(DF) >443,448,448 

443  DELF=-.75#DELF 
GO  TO  448 

445  DELF=-DEL F/2 . 

448  F=F-DELF 
DFCLD=  DF 
KTRY=KTRY+ 1 
GO  TO  425 

450  IFIKTRY-KTRYM/3 >451,452,452 

451  SF=SF/ 2 . 

452  DEL  F=S  F 

IF (KS+KK-1 >453,453,460 

453  DELYMT=YMT0UT/5. 

YMT=YMT0UT / ( 1 .-2 . * ERRM AX ) 

GO  TO  460 

455  IF(KTRY-KTRYM/3)456,456,457 

456  SYMT=S YMT/2 . 

457  DELYMT=SYMT 

IF (KS+KK-2 >458,458, 460 

458  DELF  =  (  F0UT-F)/5. 

460  KS=KS+1 

I F(KS-KSM AX >465,465,461 
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461  I F ( KFUD-KFUDM 1)462,463,463 

462  DAL=0.5*DAL 
KK=2*( KK-i )-2 
KKFUD=  1 

463  WRITE  OUTPUT  TAPE  6,464 

464  FORMAT (63H0THE  FOLLOWING  OUTPUT  FOR  THIS  LOAD  IS  OF  QUESTIONABLE  A 
1CCURACY ) 

KFUD=KFUD+ 1 

I F ( KFUD-KFUDM2 ) 500, 500,340 

465  IF ( ABSF ( I F-FOUT ) /F I-ERRM AX  1 466,466 , 47C 

466  IFIYMT 1467,470,467 

467  IF( ABSF ( I YMT-YM TOUT ) /YMT )-ERRMAX ) 500,500,470 
470  IFIKALT 1418,473,418 

473  KTRY=0 
KALT=1 

WRITE  OUTPUT  TAPE  6,474,F 

474  FORMAT! 12H  HOLDING  F  =,F12.8) 

DYMT0=YM8 

480  YMT I N=YMT 
IF(KTRY-KTRYM) 48 1,48 1,461 

481  WL=( YMB*FM(0. I-YMTJ/F 
CPL=FCP ( WL ) 

WR=(YMB*FM(  YD-YMT)  /F 
CPR=FCP!WR ) 

THETAL=FINT(FTLI»K) 

IF! IKK-2)  483,482,483 

C  IKK=2  FOR  SYMMETRICAL  LOADING,  FM<  X)*FM  IL-X) 

482  THET AR=THET AL 
GO  TO  484 

483  THETAR=FINT(FTRI,K) 

484  BCOMP=FINT(FBCI,K) 

YMTOUT=0. 

Zl= ( SUMECT  +  XNT* ! CPL-OCENT  +  FX8 ! CP  L ) ) )*(THETAL-THETAO) 
Z2=!SUMECT+XNT*!CPR-DCENT*FXB!CPR) ) ) * ( THETAR-THETAO) 
F0UT=F0+(Z1+Z2)/!CTEN+XNT*BC0MP> 

DO  485  1=1, NT 

ZZ 1= ( ECT ( I J+CPL-DCENT  +  FXBICPL) ) *  <  THETAL-THBT  AO ) 

ZZ2=(ECT(I) +CPR-DCENT  +  FXB ( CPR ) )*( THETAR-THETAO ) 

FT! I »=FTOl I)+( ZZ1+ZZ2-(F-F0)*BC0MP)/CTEN 

485  YMTOUT =YMTOUT  +  FT ! I ) *  C  ECT( I J-EC) 

DYMT=YMTOUT-YMT  IN 

490  VIN=(YMB-YMTIN)/F 

VOUT=( YMB-YMTOUT ) /FOUT 

WRITE  OUTPUT  TAPE  6 , 492 , YMT I N , V I N , YMTOUT ,DYMT , FOOT , VOUT 

492  FORMAT ( 8H  YMT IN  =,F12.8,6H  VIN  =,F1C.8,9H  YMTOUT  =,F12.8,7H  DYMT  * 
1,F12.8,7H  FOUT  =,F12.8,7H  VOUT  =,F10.8J 

IF! ABSF! ( YMTOUT- YMT IN)/YMTIN)-ERRMAX)  455,455,493 

493  IF(DYMT0*DYMT)496,455,494 
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494  IF!ABSF!DYMT0)-A8SF!DYMT) >495,499,499 

495  DELYMT=-.75*DELYMT 
GO  TO  499 

496  DELYMT=-OELYMT/2. 

499  YMT=YMT +DELYMT 
DYMTO=DYMT 
KTRY=KTRY+1 

GO  TO  480 

C  ROUTINE  FOR  DETERMINING  TENDON  STRESSES  *■  BEAM  DEFL.  1500  TO  5991 

500  IFIKK-D505,  590,505 
505  AL= ALDMBM*YMB 

WRITE  OUTPUT  TAPE  6 , 5 10, AL , YMB , F 

5  LO  FORMAT! I6H0APPLIED  LOAD  =  , F12 . 8 , 3X ,24H  MAX  APPLIED  BEND  MOM  *  ,FI 
12.8,3X,25H  TOTAL  PRELOAD  IN  BEAM  =  ,F12.8) 

DF=  F-FOUT 
DYMT=YMT-YMTOUT 
V= ( YMB-YMT ) /F 

WRITE  OUTPUT  TAPE  6 , 515 , DF , YMT ,DYMT , V 
515  FORMAT ( 7H  DIFF  = , F 1 2.8 , 3X , 21H  TENDON  BENDING  MOM  =  ,F12.8,3X,7H  DI 
IFF  =,F12.8,3X,5H  V  =  ,F12.8) 

WRITE  OUTPUT  TAPE  6 , 5 1 7, THETAL , THE T AR 
517  FORMAT (20H0LEFT  END  ROTATION  =,F12.8,21H  RIGHT  END  ROTATION  =,F12. 
18) 

WRITE  OUTPUT  TAPE  6,520 

520  FORMAT! 30H0INITIAL  TENDON  ECCENTRICITY  ,13H  TENDON  FORCE, 7X,14H  T 
1END0N  STRESS) 

DO  530  1=1, NT 
530  SIGMA!  1 1  =FT !  D/AT 

WRITE  OUTPUT  TAPE  6,540,!  EC  T!  I  ),FT(  II)  ,  SIGMA  1 1 )  ,  1  =  1  ,NT ) 

540  FORMAT! F20 . 8, 7X , 2F2C . 8 ) 

MM=MM+ 1 
TABAL! MM ) = AL 
TABTH  MM)=THETAL 
TABTR! MM)  =  THET  AR 
WRITE  OUTPUT  TAPE  6,550 

550  FORMAT! IH0,22X,6H  XD/YL,8X,11H  DEF L6CT I  ON, 5X , 1 7HCR ACK  PENETRATION) 
DO  570  I J= 1 , I JK 
XIJ=IJ 

XD=YL*XIJ/XIK 
DELT A=  PINT ! FDI , K  ) 

WXD=!YMB*FM(XD)-YMT)/F 

CPXO=FCP(WXD) 

IF(XIJ/X IK- .5 )  560,555,560 
555  TABCDt MM)=DELTA 

560  WRITE  OUTPUT  TAPE  6, 565, I J, IK, DELTA, CPXD 
565  F0RMAT(23X» 12, 2H  /, 12, 5X, 1PE15.8, 10X,0PF10.6) 

570  CONTINUE 

WRITE  OUTPUT  TAPE  6,580 
580  FORMAT! 1H  //////) 
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I F(KKFUD-1 >590,595,595 
590  F0LD=F 

YMT0LD=YMT 

YMBCLD=YM8 

595  I F ( MM— 100 ) 300t  340 »  3 AO 
END 
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CFCP  CRACK  PENETRATION  AS  A  FUNCTION  OF  W  USING  TABULATED  VALUES,  PCH 
FUNCTION  FCPIW) 

DIMENSION  FTOI 10 ) , ECT( 10 ) , FT ( 10 ) , S I GMA( 10 > , TABAL ( 100 ) , TABTL ( 100 ) , 

1  TABTRI 100), DISCI 10),NIR! 10 ), TABCP ( 200 ), TABW ( 200 >» TABCD ( 100) 

COMMON  YMB,YMT,F,F0,YL,EB,XD,WCRACK,NI,IC,I.BMD,IXS,DT0T,DCENT,DIM1 
1,DIM2,DIM3,DIM4,DIM5,DIM6,DIM7,DIM8,DIM9,DIM10, DISC, TABCP, TABW, EC 
I F ( W-WCRACK )  10,10,20 
10  FCP=0. 

GO  TO  100 
20  KL0W=1 

KH I GH=N I  ♦  1 

I F ( W-T  ABW( KHIGH ) )30,30,25 

25  FCP=DT OT 

WRITE  OUTPUT  TAPE  6,26 

26  FORMAT ( 38H0W  GREATER  THAN  WMAX,  SETTING  FCP=DTOT  ,  //  ) 

GO  TO  100 

30  KM  I D=KLOW+ ( KHIGH-KL0W)/2 
IF(W-TABWIKMID)  )  40,50,60 
40  KHIGH=KMID 
GO  TO  70 

50  FCP=TABCP(KMID) 

GO  TO  100 
60  KLCW=KMID 

70  IF(KHIGH-KLOW-l)  80,80,30 

80  SLOPE= ( TABCP ( KHIGH )-TABCP ( KLOW ))/( TABW (KHIGH J-TABW(KLOW) ) 
FCP=TABCP(  KLOW)  MW- TABW  (KLOW) >*SLOPE 
100  RETURN 
END 


C  FA  I  BEAM  CROSS  SECTIONAL  AREA  AS  A  FUNCTION  OF  CRACK  PENETRAT I  ON , PCH 
FUNCTION  FA ( CP ) 

DIMENSION  FTOI 10 ) , ECT { 10 ) , FT  I  10 ) , S I GMA ( 10 ) , TABAL 1 100) , TABTL (100), 

1  TABTR( 100) .DISC ( 10) ,NIR( 10 ), TABCP ( 200 ) , TABW (200) , TABCD (100) 

COMMON  YM8,YMT,F,F0,YL,EB»XD,WCRACK,NI,IC»IBMD,IXS,DT0T,DCENT,B,TW 
l , D , T  F , DIM5,DIM6,DIM7, DIMS, DIM9.DIM10, DISC, TABCP, TABW, EC 
IF(CP-TF)  10,10,20 
10  FA=B*( TF-CP  H-D*TW+B*TF 
GO  TO  50 

20  I F ( CP-TF-D  )  30,30,40 
30  FA=(0+TF-CP)*TW+B*TF 
GO  TO  50 

40  FA=(2.*TF+C-CP)*B 
50  IF(IXS)  70,60,70 
60  WRITE  OUTPUT  TAPE  6,65 

65  FORMAT  (33H0I  BEAM  AREA  FUNCTION  IN  PROGRAM  ) 

70  RETURN 
END 
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CFXB  I  BEAM  CENTROID  DISTANCE  FROM  CRACK  VS  CRACK  PENETRATION,  PCH 

FUNCTION  FXB(CP) 

DIMENSION  FTO( IO),ECT< 10 ) ,  FT  1 1 0 )  ,  S IGMA  UO)  ,  TAB  AL I  IOO)  ,  TABTLI  IOO  ) , 

1  TA8TRI 100), DISCI  10 ) , N IR( 10 ) , TABCP ( 200 ) , TABWI 200) , TABCD ( 100 ) 

COMMON  YMB,YMT,F,F0,YL,EB,XD,WCRACK,NI,IC,I8MD,IXS,DT0T,DCENT,B,TW 
1*D»TF»DIM5,DIM6»DIM7»DIM8,DIM9»DIM10,DI SC , TABC P , T ABW, EC 
IFICP-TF)  10,10,20 

10  T0P=0.5*B*(  TF-CP  )**2+D*TW*ITF-CP-*-0.5*D)+B*TF*(  1.5*TF-CP+D) 
FX8=T0P/FA( CP ) 

GO  TO  50 

20  IF  (CP-TF-D)  30,30,40 

30  T0P=0.5*TW*( D+TF— CP )**2+B*TF*(D+1.5*TF-CP) 

FX8=T0P/FA(CP) 

GO  TO  50 

40  FXB=(2.*TF+D-CP )/2. 

50  IF(IXS)  70,60,70 
60  WRITE  OUTPUT  TAPE  6,65 

65  FORMAT  (37H0I  BEAM  CENTROID  FUNCTION  IN  PROGRAM  ) 

70  RETURN 
END 


C  FI  I  BEAM  MOMENT  OF  INERTIA  AS  A  FUNCTION  OF  CRACK  PENETRATION,  PCH 
FUNCTION  F  I  ( CP ) 

DIMENSION  FT 0 ( 10),ECT( 10 ) , FT  1 1 0 ) , S I GMA I  10 J , TAB AL l 1 00 ) , T ABTL 1 100 ) , 

1  TABTRI 100) .DISCI 10),NIRI 10 ), TABCP I  200 ), TABWI 200 ), TABCD 1 100 ) 

COMMON  YMB,YMT ,F,F0»YL,EB,XD,WC  RAC  K,NI»IC»I BMD  »IXS»DTOT , DC ENT ,B ,TW 
1, D, TF,DIM5,DIM6,DIM7, DIMS, DIM9.DIM10, DISC, TABC P,TABW, EC 
IFICP-TF)  10,10,20 

10  F 1= I B» I TF-CP )**3+TW*D**3+B*TF**3)/12. ♦B*ITF-CP)*(FXBICP)-0.5*(TF- 
1CP) )**2+D*TW*I FXBICP )-TF+CP-0.5*D>  **2  +  B*TF*lD+1.5*TF-FXBICP)-CP) ** 
22 

GO  TO  50 

20  IFICP-TF-D)  30,30,40 

30  FI=(TW*t TF+D-CP)«*3+8*TF**3)/12.+TW*(D+TF-CP) *IFXB !CP)-0.5*(D+TF- 
1CP) )**2+B*TF*I D*1.5«TF-FXB(CP>-CP) **2 
GO  TO  50 

40  FI=(B*(2.*TF*D-CP)**3)/12. 

50  IFIIXS)  80,60,80 

60  WRITE  OUTPUT  TAPE  6,65 

65  FORMAT  (46H0I  BEAM  MOMENT  OF  INERTIA  FUNCTION  IN  PROGRAM  ) 

WRITE  OUTPUT  TAPE  6 , 70 , B , TF , D , TW 

70  FORMAT  I 16H0FL ANGE  WIDTH  =  ,F7.4,16H  FLANGE  THICK  =  ,F7.4,13H  WEB 
1DEPTH  =  , F7.4, 13H  WEB  THICK  =  ,F7.4> 

80  RETURN 
END 
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C  FI  NT  SIMPSONS  1/3  RULE  FOR  NUMERICAL  INTEGRATION 
FUNCTION  FINT(FDUM,K) 

DIMENSION  FTO( 10) ,ECT( 10 ) , FT ( 10 ) , S IGMA ( 10) , TABAL ( 100) , TABTL C 100 ) , 

1  TABTRI 100) , DISCI 10),NIR( 10) ,TABCP(200) ,TABW( 200) , TABCD ( 100 ) 

COMMON  YMB,YMT,F,FO,YL,EB,XD,WCRACK,NI, IC , IBMD , I XS ,DTOT , DCENT ,D IM1 
1,DIM2,DIM3,DIM4,0IM5,DIM6,DIM7,DIM8,DIM9,DIM10,DISC,TABCP,TABW,EC 
I1=K-1 
I 2=K-2 
XK=K 

FI NT1=0. 

FINT2=0. 

DO  10  1=1,11,2 
XI=I 

X=XI*YL/XK 

10  FINT1=FINT1+FDUM(X) 

DO  20  1=2,12,2 
X  I  =  I 

X=X I *YL/XK 

20  FINT2=FINT2+FDUM(X) 

A1=FDUM( 0.0 ) 

A2=FDUM(YL) 

FINT=( YL/( 3.*XK) )*{ Al+4. *F INT 1+2. *F INT2+A2 ) 

RETURN 

END 


CFTLI  FUNCTION  FOR  THETA  LEFT  INTEGRAL  PCH 

FUNCTION  FTLl(X) 

DIMENSION  FTO( 10),ECT( 10» ,FT( 10) , SIGMA ( 10) , TABAL ( 100) , TABTL (100), 

1  TABTR(100),DISC(10)»NIR(10)»TABCP (200) ,TABW( 200) , TABCD( 100) 

COMMON  YMB, YMT,F,FO,YL,EB,XD,WCRACK,Nl,IC, IBMD, IXS,DTOT, DCENT, DIM1 
1,DIM2»DIM3,DIM4,DIM5,DIM6,DIM7,DIM8,DIM9 ,DIM1 0 ,DI SC, TABCP, TABW, EC 
W=(YMB*FM(X)-YMT)/F 
CP=FCP(W) 

T=F#(W-EC-CP*DCENT-FXB(CP)) 

FTLI=T*( YL-X ) /( YL*EB»FI (CP ) ) 

RETURN 

END 


CFTRI  FUNCTION  FOR  THETA  RIGHT  INTEGRAL  PCH 

FUNCTION  FTRI(X) 

DIMENS  ION  FTO( 10) ,ECT( 10) , FTC  10) , SIGMA ( 10) , TAB AL ( 1 00) , TABTL ( 100 ) , 

1  TABTR( 100), DISCI  10 ) ,N IR ( 10 ) , TABCP ( 200 ) ,TABW (200) , TABCD (100) 

COMMON  YMB,YMT»F»FO»YL , EB ,XD» WCR ACK ,NI , IC , I BMD, IXS,DT0T,DCENT,DIM1 
i,DIM2,DIM3,DIM4,DIM5,DIM6,DIM7,DIM8,DIM9,DIM10, DISC, TABCP, TABW, EC 
W=(YMB»FM(X)-YMT)/F 
CP=FCP ( W ) 

T=F*(W-EC-CP+DCENT— FXB(CP)  ) 

FTRI=T*X/( YL*EB*FI(CP) ) 

RETURN 

END 


109 


PCH 


CFBCI  FUNCTION  FOR  BEAM  COMPRESSION  INTEGRAL 
FUNCTION  FBCI(X) 

DIMENSION  FTO( 10 ) , ECT ( 10 ) , FT ( 10 ) , S IGMA ( 10) » TAB AL ( 1 00) , TABTL ( 100 ) , 

1  TABTR (100 ), DISCI 10)fNIR< 10 ), TABCP ( 200 ), TABW( 200) , TABCD ( 100 ) 

COMMON  YMB,YMT,F,FO,YL,EB,XD,WCRACK,Nt , IC, I BMD , I XS ,DTOT , DCENT , D IM1 
1.DIM2, DIM3,DIM4,DIM5,D IM6,DIM7,D IM 8,0 1 M9,DI MiO , DISC, TABCP, TABW, EC 
W=(YMB*FM(X)-YMT)/F 
CP-FCP(W) 

FBCI=1./(EB*FA(CP)) 

RETURN 

END 


CFD I  FUNCTION  FOR  DEFLECTION  (DELTA)  INTEGRAL  PCH 

FUNCTION  FDIIX) 

DIMENSION  FTO( 10 ) , ECT ( 10 ) , FT( 10) , S IGMA ( 10 ) , TABAL ( 1 00) , TABTL ( 100 ) ♦ 

1  TABTR ( 100 ) t  DISCI  10 ), N IR ( 10 ), TABCP ( 200 ), TABW ( 200) , TABCD ( 100 ) 

COMMON  YMB,YMT,F,FO,YL,EB,XD,WCRACK,NI , IC, IBMD,IXS,DT0T,DCENT,DIM1 
1,DIM2,DIM3,DIM4,DIM5,DIM6,DIM7,DIM8,DIM9,DIM10,DISC,TABCP,TABW^EC 
W= ( YMB*FM ( X )— YMT ) /F 
CP=FCP ( W ) 

T=F*(W-EC— CP+DCENT-FXBICP)  ) 

G=T/(EB*FI (CP)  ) 

IF(X-XD)  10,10,20 
10  FDI=G»X*CYL-XD)/YL 
GO  TO  30 

20  FD I =G*XD* ( YL-X ) /YL 
30  RETURN 
END 


CFM  UNIFORM  BENDING  MOMENT  DISTRIBUTION  PCH 

FUNCTION  FM ( X ) 

DIMENSION  FTO(  10 ) , ECT ( 10 ) , FT< 10 ) , S I GMA < 10) ,TABAL(100) , TABTL ( 100 ) , 

1  TABTR( 100), DISCI 10),NIR( 10), TABCP! 200) ,TA8W(200) ,TABCD(100) 

COMMON  YMB, YMT ,F,FO,YL,EB,XD,WCRACK,NI, IC , I BMD , I XS ,DTOT , DCENT , D I  Ml 
i,DIM2,DIM3,DIM4,DIM5,DIM6,DIM7,DIM8,DIM9,DIM10, DISC, TABCP, TABW, EC 
RULES=*X 
FM=1.0 

IF(IBMD)  20,10,20 
10  WRITE  OUTPUT  TAPE  6,15 

15  FORMAT  ( 38H0TERM INAL  COUPLES,  C  ,  MAX  MOMENT  =  C  ) 

20  RETURN 
END 
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GLASS  I  BEAM  WITH  6  ELASTIC  TENDONS  UNDER  TERMINAL  COUPLES,  PCH 
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INCREMENT  HALVINGS  =  5  MAX  NO.  OCCURENCES  OF  QUESTIONABLE  OUTPUT 


WCRACK  =  1.041667 


TABCP 

TABW 

1 

0. 

1.04166666 

2 

0.01666667 

1.05029832 

3 

0.03333333 

1.05883005 

4 

0.05000000 

1.06725861 

5 

0.06666666 

1.07558069 

6 

0.08333333 

1.08379292 

7 

0.09999999 

1.09189187 

8 

0.11666666 

1.09987406 

9 

0.13333333 

1.10773592 

10 

0.15000000 

1.11547388 

11 

0.16666666 

1.12308426 

12 

0.18333333 

1.13056338 

13 

0.20000000 

1.13790745 

14 

0.21666666 

1.  14511262 

15 

0.23333333 

1. 15217508 

16 

0.25000000 

1. 15909091 

17 

0.26666667 

1.16585602 

18 

0.28333333 

1. 17246653 

19 

0.30000000 

1. 17891827 

20 

0.31666666 

1. 18520723 

21 

0.33333333 

1. 19132915 

22 

0.34999999 

1.19727989 

23 

0.36666666 

1.20305520 

24 

0.38333333 

1.20865084 

25 

0.40000000 

1.21406247 

26 

0.41666666 

1.21928577 

27 

0.43333333 

1.22431642 

28 

0.45000000 

1.22914998 

29 

0.46666666 

1.23378208 

30 

0.48333333 

1.23820829 

31 

0.50000000 

1.24242422 

32 

0.57500000 

1.  26084441 

33 

0.65000000 

1.27911845 

34 

0.72499999 

1.29724148 

35 

0.80000000 

1. 31520861 

36 

0.87500000 

1.  33301473 

37 

0.95000000 

1.35065430 

38 

1.02499999 

1. 36812186 

39 

1.09999999 

1. 38541155 

40 

1.17500000 

1.40251730 

41 

1.25000000 

1.41943291 

42 

1.32499999 

1.43615174 

43 

1.39999999 

1.45266713 

44 

1.47499999 

1.46897201 

45 

1.54999998 

1.48505914 

46 

1.62500000 

1.50092110 

47 

1.69999999 

1.51655009 

48 

1.77499999 

1. 53193827 

49 

1.84999999 

1.54707755 

50 

1.92499998 

1.56195973 

FIRST  DIFF 

0. 

0.00863166 

0.00853173 

0.00842856 

0.00832208 

0.00821224 

0.00809895 

0.00798219 

0.00786185 

0.00773796 

0.00761038 

0.00747912 

0.00734407 

0.00720517 

0.00706246 

0.00691582 

0.00676511 

0.00661051 

0.00645174 

0.00628896 

0.00612192 

0.00595073 

0.00577532 

0.00559564 

0.00541162 

0.00522330 

0.00503065 

0.00483356 

0.00463210 

0.00442621 

0.00421593 

0.01842019 

0.01827404 

0.01812303 

0.01796713 

0.01780611 

0.01763958 

0.01746756 

0.01728968 

0.01710576 

0.01691560 

0.01671883 

0.01651539 

0.01630488 

0.01608713 

0.01586196 

0.01562899 

0.01538818 

0.01513928 

0.01488218 
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51 

2.00000000 

1. 57657656 

0.01461683 

52 

2.07499999 

1.59091975 

0.01434319 

53 

2.14999998 

1.60498123 

0.01406148 

54 

2.22499996 

1.61875317 

0.01377194 

55 

2.29999998 

1.63222833 

0.01347516 

56 

2.37500000 

1.64540029 

0.01317196 

57 

2.44999999 

1.65826389 

0.01286361 

58 

2.52499998 

1.67081594 

0.01255205 

59 

2.59999996 

1.68305604 

0.01224010 

60 

2.67494998 

1.69498783 

0.01193179 

61 

2.75000000 

1.70662099 

0.01163316 

62 

2.82499999 

1.71797386 

0.01135287 

63 

2.89999998 

1.72907798 

0.01110412 

64 

2.97499996 

1. 73998486 

0.01090688 

65 

3.04999998 

1.75077711 

0.01079226 

66 

3.12500000 

1.76158775 

0.01081064 

67 

3.19999999 

1.77263463 

0.01104687 

68 

3.27499998 

1.78428653 

0.01165190 

69 

3.34999996 

1.79720080 

0.01291427 

70 

3  *42499998 

1.81264159 

0.01544079 

71 

3.50000000 

1.83333333 

0.02069174 

72 

3.51666665 

1. 83888887 

0*00555554 

73 

3.53333333 

1.84444444 

0.00555557 

74 

3.54999998 

1.84999998 

0.00555554 

75 

3.56666666 

1.85555555 

0.00555557 

76 

3.58333331 

1.86111109 

0.00555554 

77 

3.59999999 

1.86666666 

0.00555557 

78 

3.61666664 

1.87222220 

0.00555554 

79 

3.63333333 

1.87777777 

0.00555557 

eo 

3.64999998 

1.88333331 

0.00555554 

81 

3.66666666 

1.88888888 

0.00555557 

82 

3.68333331 

1.89444442 

0.00555554 

83 

3.69999999 

1.89999999 

0.00555557 

84 

3.71666664 

1.90555553 

0.00555554 

85 

3.73333332 

1.91111110 

0.00555557 

86 

3.75000000 

1.91666666 

0.00555556 

87 

3.76666665 

1.92222221 

0.00555556 

88 

3.78333333 

1.92777777 

0.00555556 

89 

3.79999998 

1.93333332 

0.00555556 

90 

3.81666666 

1.93888888 

0.00555556 

91 

3.83333331 

1.94444443 

0.00555556 

92 

3.84999999 

1.94999999 

0.00555556 

93 

3.86666664 

1.95555554 

0.00555556 

94 

3.88333333 

1.96111110 

0.00555556 

95 

3.89999998 

1.96666665 

0.00555556 

96 

3.91666666 

1.97222221 

0.00555556 

97 

3.93333331 

1.97777776 

0.00555556 

98 

3.94999999 

1.98333332 

0.00555556 

99 

3.96666664 

1.98888887 

0.00555556 

LOO 

3.98333332 

1.99444443 

0.00555556 

101 

4.00000000 

2.00000000 

0.00555557 
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APPLIED  LOAD  =  0.  MAX  APPLIED  BEND  MOM  =  0.  TOTAL  PRELOAD  IN  BEAM  =  12.00000000 

DIFF  =  12.00000000  TENDON  BENDING  MOM  =  0.  DIFF  =0.  V  =  0. 
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APPENDIX  B 

SELECTION  OF  A  MODELING  MATERIAL 

Because  a  meaningful  statistical  treatment  of  test  specimens 
demands  a  large  number  of  replications,  an  inexpensive  modeling 
material  was  required  for  this  program.  The  Hydrocal  gypsum 
cements  seemed  to  possess  all  the  brittle  characteristics  desired 
and  could  be  obtained  in  large  quantities  of  uniform  material  at 
a  reasonable  cost. 

To  minimize  possible  errors  in  the  measurement  of  mechanical 
strength,  the  strongest  member  of  the  gypsum  family  was  sought. 

Test  specimens  in  the  form  of  small  rectangular  beams  (l/4xl/4x4 
inches)  were  cast  using  the  Hydrocal,  Ultracal  30,  and  Hydrostone 
cements.  The  results  of  these  tests  are  shown  in  Figure  42.  The 
Hydrostone  was  selected  as  the  modeling  material  because  this  group 
(batch  A)  of  specimens  had  the  highest  average  stress. 

It  was  learned  that  the  mill  preparation  of  Hydrostone  in¬ 
volved  one  ton  lots  only.  Therefore,  to  achieve  the  desired  uni¬ 
formity  in  the  base  material  it  was  arranged  with  the  manufacturer 
to  have  the  Hydrostone  bagged  in  accordance  with  the  following 
procedure. 

(1)  Hydrostone  was  processed  in  twelve  one- ton  lots. 

(2)  Each  lot  was  kept  separate  in  the  plant. 

(3)  Two  bags  from  each  of  the  lots  were  then  selected 
and  blended  to  form  a  2400  pound  lot. 

(4)  This  process  was  repeated  to  form  ten  such  re¬ 
blended  lots. 

Following  this  procedure,  the  shipment  of  230  bags  of  Hydro¬ 
stone  were  rebagged  and  shipped  to  IITRI.  In  addition  to  the 
"A"  batch,  two  more  batches  "B"  and  "C"  were  prepared  from  two 
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different  bags  drawn  from  the  Hydrostone  supply.  The  distri¬ 
bution  of  fracture  strength  for  these  three  batches  are  shown 
in  Figure  43  .  The  small  variability  of  fracture  stress  for 
these  three  batches  indicate  that  the  precautions  taken  to  in¬ 
sure  uniformity  in  the  supply  of  Hydrostone  were  adequate. 

The  data  illustrated  in  Figures  42  and  43  were  obtained  for 
materials  mixed  at  maximum  consistency  which  was  desirable  for 
strength  purposes  but  somewhat  unsatisfactory  for  molding  pur¬ 
poses  since  the  setting  time  was  relatively  short.  For  this 
reason,  several  mixtures  of  Hydrostone  mixed  at  maximum  con¬ 
sistency  but  with  the  addition  of  minute  traces  of  sodate  re¬ 
tarder  were  prepared.  A  wet  mix  working  time  of  1-1/2  hours  was 
made  possible  by  the  following  Hydrostone-water  sodate  retarder 
ratio: 

Hydrostone  (dry):  water:  retarder  =  100:32:04.  This  formula 
was  used  uniformly  throughout  the  course  of  all  investigations 
reported  herein.  | 

Mold  preparations  and  casting  procedure  for  the  Hydrostone 
specimens .  -  Various  component  shapes  were  required  throughout 

the  experimental  effort.  These  shapes  ranged  from  simple  beams 
of  square  cross  section  to  shapes  of  a  somewhat  more  complicated 
geometry  involving  surfaces  curved  in  two  directions.  The  simple 
shapes  such  as  the  2.5  inch  x  2.5  inch  x  40  inch  beam  and  the  15  inch 
diameter  by  5/8-inch  thick  plate  were  cast  in  one  piece  molds.  The 
mold  for  the  beams  was  constructed  as  a  battery  of  five  units.  See 
Figures  44  and  45  .  The  complicated  mold  configurations  as  re¬ 
presented  by  the  "backbone"  and  the  ogive  shell  were  obtained  by 
first  making  master  molds  using  wood  and  aluminum  materials.  The 
master  molds  were  then  used  to  prepare  room- temperature  vulcani¬ 
zing  rubber  ship  molds.  These  molds  were  in  turn  used  as  the  pro¬ 
duction  mold  in  casting  the  required  quantity  of  the  desired 
shape.  See  Figures  46  ,  47  ,  and  41  . 
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Figure  43  Distribution  Curves  of  Fracture  Stresses  for  Three  Batches 
of  Hydrostone  Gypsum  Cement 


Figure  44  Example  of  Battery  Mold  for  Casting  Beams 
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Figure  46  Example  of  Cylinder  and  Backbone  Specimens 


The  following  casting  process  was  used  to  produce  consistent 
Hydrostone  plaster  test  specimens;  Hydrostone,  water,  and  retarder 
were  weighed  out  in  the  normal  consistency  proportions  of  100:32:04. 
The  gypsum  was  slowly  sifted  into  the  water  and  allowed  to 
soak  for  about  one  minute.  The  mix  was  then  stirred  by  hand  tor 
three  to  five  minutes  followed  by  mechanical  mixing  for  three  to 
four  minutes.  The  mechanical  mixing  coalesced  the  smaller  bubbles 
into  larger  bubbles  which  quickly  came  to  the  surface  of  the  mix 
and  dissipated  into  the  air.  The  resulting  castings  were  reasonably 
bubble  free  and  fairly  consistent  in  appearance  from  batch  to  batch. 

Curing  of  Hydrostone  specimens.  -  In  order  to  insure  uniformity 
among  specimens  with  regard  to  strength,  each  specimen  was  cured 
by  driving  off  the  excess  moisture  in  a  drying  oven.  It  was  deter¬ 
mined  that  in  the  case  of  the  large  volume  beam  specimens,  ten 
hours  of  exposure  at  104°F  was  sufficient  to  bring  the  specimens 
to  constant  weight.  As  an  added  margin,  no  specimen  tested  was  sub¬ 
jected  to  less  than  24  hours  of  drying  time. 
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